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ABSTRACT 

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long 
narrow streams of merging galaxies and a smoother component of cold gas. We post-process 
seven high-resolution simulated galaxies with radiative transfer to study the absorption char- 
acteristics of the gas in galaxies and streams, in comparison with the statistics of observed 
absorption-line systems. We find that much of the stream gas is ionized by UV radiation from 
background and local stellar sources, but still optically thick (iVni > 10^^ cm^^) so that the 
streams appear as Lyman-limit systems (LLSs). At z > 3, the fraction of neutral gas in streams 
becomes non-negligible, giving rise to damped Lyman-a (DLA) absorbers as well. The gas 
in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within 
one (two) virial radii, the covering factor of optically thick gas is < 25% (10%) for LLSs 
and < 5% (1%) for DLAs, slowly declining with time following the universal expansion. 
Nevertheless, galaxies and their cold streams in the studied mass range. Mm — 10^° — 10^^ 
Mq, account for > 30% of the observed absorbers in the foreground of quasars, the rest 
possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in 
the streams is ~ 1% solar, much lower than in the galaxies. The simulated galaxies repro- 
duce the Lya-absorption equivalent widths observed around Lyman-break galaxies, but they 
severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise 
from outflows. We conclude that the observed metal-poor LLSs are likely detections of the 
predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations 
with more massive outflows in conjunction with the inflows may enable a clearer distinction 
between the signatures of the various gas modes. 

Key words: Radiative transfer; galaxies:formation; galaxiesievolution; galaxies:high- 
redshift; quasars:absorption lines; intergalactic medium. 



1 INTRODUCTION 

Cold neutral hydrogen that fuels star formation is an important 
driver of galaxy evolution at all times. Recent progress in ob- 
servations and simulations has remarkably improved our view of 
gas in the early Universe. Molecular gas in emission is detected 
in normal star for ming galaxies to z ~ 2.3 l lTacconi et alj|201ol : 
lOaddi et"aill201(l) . Large imaging and spectroscopic surveys have 
increased by orders of magnitude the samples of quasars and galaxy 
pairs useful for probing in absorption the interstellar medium 
(ISM), the circumga lactic medium (CGM) and the intergalactic 
medium (IGM) (e.g. lAdelberger etai]|2003l : iHennawi et alll2006l : 
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20ld:ISteidel et 31.1201 0*). Similarly, the increased resolution of nu- 
merical simulations has made it possible to repr oduce the morphol- 
ogy of galaxies especially at high redsh ift (e.g. lAgertz et al.ll2009l: 
iDekel et alj2009a| - |Ceverino et al.ll201(]t) . 

Significant progress is being made on the basic issue of 
how gas is accreted from the IGM into galaxies. A consis- 
tent theoretical picture i s emerging from both analytic theory 
( iBimboim & Dekell l2003l: [Dekel & Bimboim 2006) and numeri- 
cal simulations (Keres et alj |2005l : iGcvirk et alj|200^ : iKeres et all 
[2009; Dekel et al. 2009b). High-redshift galaxies are predicted to 
acquire the bulk of their gas mass at temperatures of a few lO* K 
via long narrow streams that consist of a smooth component and 
small merging galaxies. However, indisputable observational evi- 
dence for this predicted dominant mode of gas input is still lacking. 



2 Fumagalli et al. 



Lya emission due to cooling radiation has been proposed 
as a signature of cold gas accretion that powers Lya blobs 
(e.gJ Haiman et al."2000"; Fardal et al."200l'; Furlanetto et al.'ioO?; 
iDiikstra & Lo eb 2009; Goerdt et al. 2010. ; Matsuda et al. 2010), 
but the resonant nature of the Lya line and the different 
physical mechanisms that i nfluence the total Lya luminosity 
jpaucher-Giguere et al.l l2010l Kasen et al., in prep.) prevent us 
from unambiguously connecting the observed fluxes with the pres- 
ence of cold streams. Observations of gas in absorption against 
background sources are, in principle, ideal to probe the COM to 
uncover infalling gas. However, it has been suggested that the 
low covering factor of the cold filaments reduces the probability 
of detecting t hem in absorption jpaucher-Giguere & Kereslboill ; 
IStewart et al.l 2010i), espe ciallv in the presence of large scale out- 
flows dSteidel et al.ll2010l) . 

The subject of the present paper is the problem of detecting 
streams of cold gas in absorption against background sources. We 
formulate detailed theoretical predictions of the covering factor, the 
cross section, the column density distribution, the kinematics, and 
the metallicity of the neutral gas in seven galaxies, simulated with 
a high-resolution AMR code in a cosmological context and ana- 
lyzed in the redshift range z = 1.4 — 4. To facilitate the compari- 
son between our theoretical predictions and observations, we adopt 
the same methodology used in observational studies of absorption 
line systems (ALSs) in background quasars and galaxy pairs. This 
work focus on the visibility of cold stream and by comparing and 
contrasting these predictions with observations, we aim to find the 
observables that are best suited to probe this mode of accretion. 
A companion paper (Kasen et al., in prep.) will focus on the Lya 
emission properties using the same simulations. 

The paper consists of three parts. In the first two sections, we 
briefly discuss the numerical simulations (Section |2j and the ra- 
diative transfer calculation used to estimate the hydrogen ioniza- 
tion state (Section |3]l. This part is complemented with Appendix 
lAl where we provide detailed comparisons of the neutral hydro- 
gen masses and column density distributions from three different 
radiative transfer models. In the second part, we quantify the neu- 
tral hydrogen cross section and covering factor as functions of the 
gas column density and redshift (Section 0. In the third part, we 
compare more directly model predictions with observations, focus- 
ing in particular on the column density distribution and incidence 
of absorption line systems (Section |5} as well as their metallicity 
and kinematics (Section|6}. A summary and conclusions follow in 
Section|7] 

Throughout this wor k, we assume the W MAP5 values of the 
cosmological parameters dKomatsu et alj2009l) . consistent with the 
cosmology adopted in the simulations (f2m = 0.27, Q,h ~ 0.73 
and //o = 70 km s"^ Mpc"^). 



2 THE SIMULATIONS 
2.1 Simulation method 

A more detailed discussion of the simulations used in this work and 
comparison s with properties obse rved in high-redshift galaxies can 
be found in lCeverino et al.l ( l20l tf). Here, we briefly summarize the 
numerical procedures adopted and we review some of the uncer- 
tainties associated to these simulated galaxies that are relevant to 
this analysis. 

The simulations have been run with the adaptive mash refine- 
ment (AMR) hydro-gravitational code ART jKravtsov et al.lll997l ; 



lKravtsovll2003h . The code incorporates many of the relevant physi- 
cal processes for galaxy formation, including gas cooling and pho- 
toionizatio n heating, star f ormation, metal enrichment and stellar 
feedback ( Ceverino & Klvp in 2009; Ceverino et al. 2 010) . Cooling 
rates were tabulated for a given gas density, temperat ure, metallic- 
ity, a nd UV background (UVB) based on CLOUDY jperland et"^ 
1 19981) and assuming cooling at the center of a one kpc thick cloud, 
illuminated by a uniform radiation field. Photoheating is self- 
consistently modeled together with radiative co oling, through the 
same CLOUDY tables, assuming a uniform UVB ( IHaardt & Madaul 
1 19961) . Self-shielding of dense, galactic neutral hydrogen from the 
UVB is approximated by suppressing the UVB intensity to the pre- 
reionization value (5.9 x 10~^^ erg s~^ cm^^ Hz^^ sr^^) above 
gas densities nu —0.1 cm~"^. This threshold is consistent with the 
results from radiative transfer calculations (see Appendix [At. 

Stars form according to a stochastic model that is roughly con- 
sistent with the KennicutI ( 1998) law in cells where the gas temper- 
ature is below lO^K and the density is above the threshold nn ~ 1 
cm~^. The code implements a feedback model, in which the com- 
bined energy from stellar winds and supernova explosions is re- 
leased at a constant heating rate over 40 Myr, the typical age of the 
lightest star that explodes as a type II supernova (SNII). Pnergy 
injection by type la supemovae (SNIa) is also included. The heat- 
ing rate from SNIa assumes an exponentially declining rate with a 
maximum at 1 Gyr. 

The ISM is enriched by metals from SNII and SNIa. Metals 
are assumed to be released from each stellar partic le by SN II at a 
consta nt rate for 40 Myr since its birth, assuming a lMiller & Scald 
(Il979h initial mass function (IMP) and matching the yields of 
IWooslev & Weaved ^995). The metal ejection by SNIa assumes 
the same exponentially declining rate. The code treats metal ad- 
vection self-consistently and metals can diffuse and pollute the 
medium around the central disks. In these simulations, the dark 
matter particle mass is 5.5 x 10® Mq, the minimum star particle 
mass is 10* Mq, the smallest cell physical size is always between 
35 and 70 pc. 

Pach dark matter halo has been selected from a cosmological 
dark-matter-only simulation and resimulated with the addition of 
baryons, using a zoom-in technique that follows the evolution of 
the selected lagrangian regi on of ~ 1 comoving Mpc ^ (for details 
on the initial conditions see lCeverino & KlvpirjhoogI) . All the sim- 
ulated galaxies are centrals to their halos and the analysis has been 
performed on a box of four times the virial radius (i?vir), centered 
at all times on the main halo. Note that the selected galaxies do 
not reside in a sizable group or cluster at z = because of the 
difficulties in simulating companion galaxies. 

Although these simulations reproduce some of the main prop- 
erties observed in high-redsh ift galaxies, like the T ully-Pisher and 
the mass-SPR relations (see ICeverino et al. I I2OIO '), the use of ap- 
proximated prescriptions in lieu of the complex interplay of physi- 
cal processes that operate on scales smaller than our resolution lim- 
its the predictive power of these models and of simulations in gen- 
eral. Particularly relevant for our work is the uncertainty associated 
with stellar feedback, imp lemented according to the prescription of 
ICeverino & Klvpinl ( 120091 ). In these simulations, outflows eject hot 
gas and metals at hi gh velocity outside the ce ntral disks into the 
IGM and CGM (see ICeverino & KlYpid l2009'). The total outflow 
mass flux reaches, in some cases, about one-third of the inflow flux. 
However, th ese winds may not be as strong as some observations 
suggest (e.g. lSteidel et all2010l) and AGN feedback is missing. As 
a consequence, a fraction of the gas inside the disks might actu- 
ally be ejected by stronger outflows, altering the covering fraction 
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Figure 1. Evolution of the galaxy properties within the virial radius. Each galaxy is identified by a different line connecting different symbols, a) Virial radii, b) 
Dark ma tter mass, where the dotte d line marks the characteristic halo mass, and the horizontal dashed line is the critical mass above which a stab le shock can 
develop jPekel & Bimboin]|2006l) . c) Star formation rate, where the shaded area represents the typical SFR of Lyman break galaxies at z ~ 2 jReddv et alj 
I2OI0I) to highlight the general overlap, d) Stellar masses, e) Hydrogen masses, f) Gas fraction, defined as the ratio of gas mass to total baryon mass. 



of neutral gas and the metall icity distribution in the CGM (see e.g. 
iFaucher-Giguere et al.ll201lh . Unless the low density winds dras- 
tically perturb the shape and kinematics of the gas that inflows in 
dense narrow streams, these models capture the properties of the 
cold streams, but caution is advised in generalizing our results. 



2.2 Sample properties 

For each galaxy, the total dark matter and instantaneous stellar mass 
are computed summing all the particles within the virial radius. 
Similarly, we obtain the total star formation rate (SFR) within the 
virial radius by averaging the masses of all the stars formed in the 
last 60 Myr. This time interval is arbitrarily chosen to be large 
enough to avoid fluctuations due to the stochastic star formation 
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rate implemented in the simulation, but the average SFR does not 
depend on this assumption. The total hydrogen masses are com- 
puted from the hydrogen number density in each cell of the AMR 
grid as mn = irip tih , with rrip the proton mass and / the cell 
size. No correction for heavier elements, including helium, is ap- 
plied throughout this work, except when we quote a gas fraction 
where we multiply the hydrogen mass by 1 .38 to correct for heavier 
elements. Finally, we derive virial masses by combining together 
the baryonic and non-baiyonic masses. 

The sample is composed of seven massive galaxies that ex- 
hibit dark matter and baryonic properties roughly consistent with 
the observed population of Lyman break galaxies (LBGs). Galaxy 
properties as a function of tim^^are displayed in Figure[T|and listed 
in Table lBll and Table lB2l in Appendix|B] Panel a) and b) show the 
redshift evolution of the virial radius and dark matter mass. These 
galaxies span roughly more than a decade in mass. As seen in panel 
b), they are all well above the characteristic Press-Schechter halo 
mass (M») in the relevant mass range, more so at higher redshifts. 
Galaxy MWl, and to some extent MW3, exhibit a faster overall 
growth rate than the other galaxies. The interval in mass and red- 
shifts analyzed here is typical for cold gas accretion; only at later 
times {z < 2) halos approach the critical mass (horizontal dashed 
line) above which a stable shock develops and cold accretion b e- 
comes less efficient jPekel & BimboinJl2006l ; lKeres et all2009l) . 

The evolution of the SFR within the virial radius is presented 
in panel c). While there is a tendency to have higher star formation 
rates in higher mass halos, the stochastic nature of the simulated 
SFR and different accretion histories produce a less ordered star 
formation history (SFH). The star formation in the sample brack- 
ets the typical va lues found in LBGs (shaded area for z ~ 2; 
iRedd v et al.' 2010), with MW8 and MW9 usually below the com- 
monly observed range. The evolution of the stellar mass is shown 
in panel d) and follows the general halo growth. Finally, the gas 
properties are displayed in panels e) and f) where we show the time 
evolution of the hydrogen mass and of the gas fraction, here defined 
as the ratio of gas mass to total baryonic mass. The simulated galax- 
ies have a gas fraction com parable to the estimates from observe d 
molecular gas at 2 ~ 1 - 2 ffacconi et alj2010l : lDaddi et alj201oh . 
The baryonic fraction inside the virial radius is at all times below 
the universal value /b = 0.165, approaching 0.25 — 0.3/b when 
only gas and stars are considered inside the central disks. MW2 is 
an exception, having a particularly high baryon fraction at all red- 
shifts, probably reflec ting an extraordinary high gas accretion rate 
( ICeverinoetal.l2010l) . 



2.3 Satellite identification 

Being part of a larger cosmological simulation, each box contains 
also satellites whose positions and sizes are used during the anal- 
ysis to disentangle gas that resides in bound structures from gas 
in the streams. For t his reason, we identify the satellites using th e 
AMIGA halo finder jGill et al.ll2004l : iKnollmann & Knebdl2009l) . 
This code locates the centers of the halos and iteratively computes 
their radii so that they define an overdensity Avir over the back- 
ground density ph- For small child halos embedded within larger 
parent halos, this radius is defined instead as the distance at which 
the density profile reaches a minimum as one moves away from the 
halo center We adopt the internally computed Avir from AMIGA as 



^ Galaxies MW2, MW4, and MW5 were not available at low redshifts 
when this sample was assembled. 



a function of redshift. We consider particles with velocities greater 
than the escape velocity to be unbound and we select only satellites 
with more than 182 particles. This sets a lower limit to the satellite 
halo mass of ~ 10* Mq. 

While AMIGA outputs several parameters for each of the iden- 
tified halos, in the following we compute masses and SFRs con- 
sidering all the particles that are enclosed in the defining radius, 
purely based on positions. Although not identical to the mass com- 
puted internally by AMIGA, these masses are well correlated with 
the ones generated by the code. This simple definition allows us 
to compute self consistently the Hi masses and SFRs. Due to this 
approximation, our satellite masses can be uncertain up to a factor 
of two. Child halos contribute to the parent halo masses. 



3 ATOMIC HYDROGEN NEUTRAL FRACTION 

3.1 The radiative transfer post processing 

A reliable estimate of the ionization state of the gas is essen- 
tial for any study of absorption line systems. Ideally, the ion- 
ization state of the gas should be coupled to the hydrodynamic 
calculations, but due to the numerical complexity of this prob- 
lem, codes tha t solve the (approxim a ted) radiative transfer (RT) 
equation (e.g. iGnedin & Abell 1200 ll : lAbel & Wandeli l2002h to- 
gether with the hydrod ynamic equations has been developed and 
used only recently (e.g. Razoumov et al ]|2006l : |Petkova & Springell 



|2009|; IWise& Abelll2oT ll) and they remain expensive for high res- 
olution simulations. Currently, RT post-processing is the practical 
way to compute the ionization state of the gas in our high-redshift 
simulations. The post processing approximation is justified by the 
fact that the radiation-diffusion timescale is short compared to the 
evolution timescale. 

The neutral fraction xhi of the atomic hydrogen in each AMR 
cell is determined using a Monte Carlo RT code that includes both 
coUisional ionization and photoionization due to Lyman contin- 
uum radiation from both the cosmological UVB and local stellar 
sources, together with absorption and scattering by dust grains. 
This calculation improves upon results based on simple ray-tracing 
schemes. In fact, Monte Carlo transport does not introduce pre- 
ferred directions in the photon path and takes into account the scat- 
tering of ionizing photons. Thus, a diffuse component of the UV ra- 
diation field is naturally included (about 32% of the time, photons 
scatter off of hydrogen atoms rather than being absorbed). Further, 
radiation from stars and dust absorption is modeled locally for each 
stellar particle in order to reproduce the anisotropy of the radiation 
field within the central disks and satellites galaxies. Additional de- 
tails on the numerical procedures and validating tests can be found 
in a companion paper (Kasen et al., in prep.). Here, we briefly sum- 
marize the models adopted for the source and sink of ionizing radi- 
ation. 



3.1.1 CoUisional ionization 

At temperatures above lO** K, hydrogen atoms are ionized by elec- 
tron collisions. Neglecting cosmic rays that can affect the neutral 
fraction even at lower temperatures, electron collisions are the only 
factor that regulates the ionization state of the gas in the absence 
of external radiation fields or in regions that are optically thick to 
UV radiation. In our simplest model, we assume coUisional ion- 
ization at equilibrium (CIE) under the si mplifying assumption of a 
dust free gas without external radiation dOnat & Stemberg]|2007h . 
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Figure 2. Ionizing radiation and dust in tlie central disk of MW3 (face on). Left panel: surface density of ionizing radiation smoothed witli a Gaussian kernel 
of 240 pc, showing giant clumps with enhanced star formation. Right panel: integrated optical depth of dust at 912A along the entire depth of the box (282 
kpc). Star formation occurs primarily in overdense and enriched regions where the dust optical depth is high. 



Under this approximation, gas below T — IQ K is fully neutral, 
while gas above T > 2 x 10* becomes highly ionized. In models 
that include also UV radiation, collisional ionization is coupled to 
photoionization. Since our RT calculation does not include photo- 
heating, the final temperature and the resulting ionization by colli- 
sions are likely to be underestimated in partially shielded regions. 
This problem is alleviated by the fact that the hydrodynamic code 
already includes a treatment for photo-heating Q 



3.1.2 UV radiation from cosmological background 

The UV radiation from unresolved background sources is the next 
ingredient we add to the RT calculation. Since the UVB spec- 
trum does not vary much with frequency over the relevant inter- 
val for hydrogen ionization, we assume a constant mean intensity 
J912 = 5 X 10"^^ erg s"^ cm~^ Hz"^ sr"^ ^aardt & Madau 
over the energy range 13.6 eV — 54.4 eV, corresponding to 
the H I ionization threshold and the cutoff at the He 11 threshold, re- 
spectively. Further, the modest variation as a function of redshift in 
the i nterval g ^ 2 — 4, wi thin the un certainties of the UV B models 
(e.g. lFaucher-Giguere et a l. 2008; Pall' Aglio etd]|2009l) . justifies 
the use of a fixed value with time. Our analysis is sensitive to the 
specific assumptions of the UVB, but it will become evident from 
the discussion that a decrease in J912 from 2 = 3 to 2 = 4 would 
not significantly affect (and even reinforce) our results. 



3.1.3 UV radiation from local sources 

An additional and important source of local ionizing radiation is 
the Lyman continuum from newly born stars. In our model, we es- 
timate the rate of ionizing photons {Qua) emitted by a young star 
using stellar population synthesis models. Qjjo is a function of the 
metallicity, IMF, SFH, and age of the stellar populations. Among 
those, the most dramatic dependence is on the SFH and age, since 
the ionizing radiation from a single burst can vary up to several or- 
ders of magnitude in a short time interval (1 — 100 Myr), when the 
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most massive stars that contribute to the Lyman continuum lumi- 
nosity leave the main sequence. 

As a first approximation, we estimate the total number of 
ionizing photons from a constant mean SFR from the entire sim- 
ulated box (i.e. from both the central galaxy and satellites) us- 
ing the proportionality (5jjo(s~^) — a SFR{MQ/yj:~^). We 
compute the normaliz ation log a = 53.168 using STARBU RTS99 
( lLeithereretai1ll999h . with a Kroupa IMF jXroupjlloOll) above 
0.08 Mp) , su bsolar metallicity Z = 0.004, and stellar models from 
iLeieune & Schaerer (2001). Stellar particles have different degrees 
of enrichment and the mean metallicity ranges from Z = 0.004 to 
Z = 0.013 in our sample, but the variation of Q^o with metallic- 
ity is at most a factor of two across the range Z = 0.001 — 0.04 
(log a = 53.227 - 52.935), and even less (~ 0.1 dex) with the 
choice of the available stellar models. We redistribute the total Qpjo 
among stellar particles which are younger than 20 Myr, in a fraction 
that is proportional to the stellar mass Mi^stm 



.HO 



M,,,tarft < 20 Myr) 
Y:,M,,,t^.r{t < 20 Myr) 



(1) 



Beside the intrinsic uncertainties in the stellar models and 
metallicity, other reasons make the photoionization rate from stars 
difficult to estimate. For example, our finite resolution affects the 
ability to properly resolve the clumpy ISM at scales below 35 — 70 
pc. This enhances the leakage of Lyman continuum photons due to 
the absence of high density structures around stellar particles, but at 
the same time the escape fraction of ionizin g photons is decreased 
by the lack of porosity jWood et"al]|2005l) . A proper estimate of 
the error associated with the leakage of UV photons is difficult, 
but most of the ionizing radiation from stellar sources is absorbed 
in the proximity of the emitting regions at high hydrogen volume 
density. In these models, the escape fraction at the virial radius is 
between 5 — 10 %, varying with redshift and from galaxy to galaxy. 
As consequence, only a small fraction of the escaping radiation ion- 
izes gas near the central disks with A'^hi ~ 10^* — 10^" cm~^ (see 
next section). Therefore, the associated uncertainties on the neutral 
gas and the cumulative cross section of optically thick hydrogen 
beyond 7?vir with TVhi > 10^^ cm~^ are expected to be minor. 

The left panel of Figure |2] shows the ionizing radiation rate 
surface density (convolved with a Gaussian kernel of 240 pc only 
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Figure 3. Hydrogen column density for MW3 at z = 2.3. Top left: A'^h- Top right: Afjj i from the CIE model. Bottom left: Nji i from the UVB model. Bottom 
right: A'^h I from the STAR model. Most of the gas that resides in the streams is ionized by electron collisions and the UVB, while photons from newly bom 
stars affect the high column density inside the central and satellite galaxies and their immediate sun'oundings. 



for visualization puiposes), for the central disk of MW3 at z — 
2.3 and in a nearly face-on view. Most of the ionizing radiation is 
produced in clumps and in spiral density wakes. While the SFR in 
the entire box approaches 50 Mq yT~^, the central disk accounts for 
only half of the total ionizing radiation (with SFR ~ 27 Mq yr~^). 
About 50% of the remaining photons come from a large satellite 
which is forming stars at ~ 10 M0 yr^^. This example highlights 
the fact that photons from young stars that escape the surrounding 
halos can introduce low level anisotropics in the UV radiation field. 



3.1.4 Dust opacity 

Dust is an important sink of Lyman continuum photons, particu- 
larly relevant for radiation from local sources. Using the metallicity 
in the gas phase from SN II and SN la, the dust volume density is 
given by 



Pd = /d M mp (ZsN la + ^SN ii)(nHi + 0.01 nHii) 



(2) 



where fd = 0.4 is an estimate of the fraction of metals locked 
in dust i lDweklll998l) and /X = 1.245 is the mean particle weight, 
including helium. To mimic grain destruction, we suppress dust for- 
mation in ionized regions to only 1 % of what is found in the neutral 



pha se. This fraction is th e most uncertain quantity in the dust model 
(see lLaursen et al .1200 9', section 7). 

From the dust density we compute the dust optical depth at 
912A, i.e. the hydrogen ionization potential, in each AMR cell as 
r = a I, with a — Ud (cts + fa) = k pd (1 — and I the 

cell size. Here, as and (Ta are the scattering and absorption cross 
sections, A is the albedo, and k is the frequency dependent dust 
absorptive opacity. For the above quantiti es, we assume k — 9.37 x 
10'' cm^ g"^ and A = (Js/ad = 0.248 ( lDrainell2003h . During the 
RT calculations, at the relevant UV wavelengths, a linear fit to the 
iLi & Draind ( I2OOII) data is used: k{\) = 9.25 x lO'' - 91.25 x 
(912A - A) and A(A) = 0.24 + 0.00028(912A - A) . A map of the 
dust optical depth for MW3 at 2 = 2.3 in a nearly face-on view is 
in the right panel of Figure|2] Here the projected dust optical depth 
is computed along a path of 282 kpc, the size of the entire box. 



3.2 Results of the radiative transfer calculation 

For each galaxy, we run three different RT models, gradually in- 
cluding additional physical processes. In the first model (hereafter 
CIE model), we derive the neutral fraction assuming CIE, without 
any source of radiation. In the second calculation (UVB model), 
we include the UVB together with dust and collisional ionization. 
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Finally, in our third model (STAR model) we add ionizing radia- 
tion from local sources to the UVB model. Figure |3] presents an 
output from these calculations. In the top left panel, we show the 
projected A^h column density for MW3 at 2 = 2.3. High column 
density gas is accreting onto the central galaxy through large ra- 
dial streams, with gas overdensities associated with satellites (two 
incoming galaxies along the streams and two closer in, near the 
central disk). In the other three panels, we display the A^h i column 
density from the different RT models. 

This figure captures the basic differences that arise from the 
different physical processes included in the RT calculation. Part of 
the gas within the filaments has a temperature above ~ 10* K, and 
coUisional ionization alone (top right panel) lowers the neutral col- 
umn density by more than one order of magnitude. A comparison 
of the H I map for the CLE and UVB models (bottom left) clearly 
shows that the CIE approximation largely overestimates the neutral 
fraction and that photoionization from the UV background affects 
most of the low density gas in the streams. Indeed, the filaments 
are highly ionized, with patches of self-shielded neutral gas that 
surround the main galaxy and the satellites. Cold streams are not 
entirely self-shielded. Finally, the inclusion of local sources mostly 
affects the high column density gas (where stars form) and their im- 
mediate surroundings, where the column density is caused to drop 
below A'^Hi ~ 10^" cm~^. The low escape fraction from the galaxy 
disks (below 10% at the virial radius) implies a minor effect on 
the gas in the streams beyond i?vir without any appreciable dif- 
ference compared to the UVB model for column densities below 
A^Hi ~ 10^* cm^^ 

A more quantitative comparison between the different models 
is presented in Appendix lAl There, we discuss the typical volume 
density for self-shielding, and the effect of local sources on the 
column density and mass of neutral hydrogen. We also provide a 
crude fitting formula to the UVB model useful to improve the CIE 
approximation. 



4 THE PROJECTED AREA OF NEUTRAL GAS 

In the previous section, we have shown that the streams are only 
partially shielded and that both the UVB and local sources are nec- 
essary ingredients to study the neutral and ionized gas distribution 
across a large range of column density. In this section, we quan- 
tify the cross section of neutral gas in and around massive galaxies 
and the evolution of the covering factor. In the remaining of our 
analysis, we consider the STAR model as our fiducial RT calcula- 
tion, providing from time-to-time a comparison with results from 
the UVB model. 



4.1 Neutral hydrogen cross section 

We derive the cross section that these simulations offers in projec- 
tion as a function of the H I column density by collapsing the AMR 
3D grid onto a 2D regular grid with cell size equal to the size of the 
finest level in the AMR structure. We measure in each galaxy the 
covered area in A^h i bins of 0.2 dex, averaged in three orthogonal 
directions along the axes of the AMR grid. An example for MW3 
at2 = 2.3 (STAR model!) is in Figure|4] For relative comparisons, 
we introduce in the right hand side axis the covering factor, defined 



^ We highlight how the three RT calculations affect the cross section at 
various column densities in Figure IaTI 
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Figure 4. Cumulative cross section for MW3 at 2 ~ 2.3 (STAR model) 
averaged along three orthogonal directions. The different curees are for a 
cylinder of radius 2i?vir(solid line), a cylinder of radius 0.2i?vir (dotted 
line) and for the streams only within 2i?,vir (dash dotted line). The hatched 
region between dashed lines highlights the dispersion along the three pro- 
jection axis for a cyhnder of radius 1 -Rvir ■ On the right axis, the physi- 
cal area is translated into a fractional covering factor, normalized to the 
area within 2/?vir- Colors represent intervals of column density for DLAs, 
SLLSs, LLSs and MFP gas. Gas within galaxies contributes to half of the 
DLA cross section, while LLSs and MFP gas are mainly associated with 
cold streams. Only ~ 10% of the projected area is covered by optically 
thick gas and ~ 1% of the cross section is covered by primarily neutral gas 
at A^HI > 10^" cm-2). 



as the cross section normalized to the area inside a circle of radius 
2_Rvir. The solid line represents the cross section measured in the 
entire box, while the dot-dashed line highlights the contribution of 
the cold streams alone, obtained by masking all the gas cells in- 
side one quarter of the virial radius of the main halos and satellites. 
Further, to isolate the contribution of the central galaxy, we display 
the cross section within two concentric cylinders of sizes 0.2_Rvir 
(dotted line) and _Rvir (dashed lines). In the latter case, we display 
the minimum and maximum cross section among the three projec- 
tion axis with a hatched region. Obviously, different orientations 
have different cross section distributions, but the dispersion along 
the mean appears not to exceed the amplitude of the features that 
are visible in the distribution. The radius 0.25-Rvir is well suited 
to separate streams and galaxies and it encompasses the interphase 
region where the flowing gas is recoimecting with the central disk. 
However, the fraction of neutral gas in the column density interval 



10^ 



10 cm that is associated with the streams depends on 



the location of this boundary. 

We mark in this figure four A*'h i intervals that define the ion- 
ization state of the atomic hydrogen and that are associated with 
different classifications of observed ALSs. This is done to high- 
light at what intervals of column densities cold streams are pre- 
dicted to dominate the neutral gas cross section. Specifically, at 
A^'hi ^ 2 X 10^° cm~^, the hydrogen is neutral, highly optically 
thick (r > 1000), and gives rise to damped Lyman-a absorbers 
(DLAs). The interval 10^^ cm'^ A^hi < 2 x 10^° cm"^ marks 
the transition between fully neutral to ionized hydrogen and defines 
the super Lyman limit systems (SLLSs) or sub-damped Lyman-a 
absorbers. The interval 1.6 x 10^^ cm"^ < A^hi < 10^^ cm"^ 
defines Lyman limit systems (LLSs) which are highly ionized, but 
retain enough neutral hydrogen to be optically thick (1 < r < 60). 
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Figure 5. Average DLA cross section within -Rvir as a function of halo 
mass for the central galaxies only (blue squares) and central galaxies plus 
satellites (black circles). All the galaxies in our sample, at all sampled red- 
shifts, are included. The vertical error bar indicates the standard deviation 
about the mean, while the horizontal error bar indicates the interval of halo 
mass used in the average. Also shown, individual galaxies at z ^ 2.8 (green 
downward triangle) and z ~ 1.9 (red upward triangle). A linear regression 
(dashed line) shows that the total cross section is nearly proportional to 
the halo mass below 4 X 10^ Mq, while it is roughly constant at higher 
masses. A shallower dependence is found over the entire available mass 
range (dotted line), consistent with DLAs originating above a fixed suiface 
density threshold for a self-similar gas density distribution. Central galaxies 
typically account for less than 50% of the DLA mean cross section. 



Finally, below A^hi ~ 1-6 x 10 cm~ , the hydrogen becomes 
optically thin (t < 1) and it is highly ionized. However, this gas 
is estimated to contribute to the mean free path (MFP) of Lyman 
photons in the Universe at z > 3 |Prochaska et al. 2009, 2010). In 
this work, we consider the interval 3.2 x 10^^ cm"'^ ^ A'^hi < 
1.6 X 10^^ cm~^, which we dub as "mean free path" gas. 

The cross section distribution within iZvir steeply rises from 
the high end towards A^hi ~ 10^*^ cm~^ and then flattens between 
iVm ~ 10'* - 10^° cm where gas is ionized. At lower column 
density there is a mild increase. Gas within IRvii exhibits a sim- 
ilar shape in its cross section, but the separation between the two 
distributions increase systematically moving towards low column 
densities. From the example provided in Figure |4] we find that at 
2-Rvir, ~ 10% of the projected area is covered by optically thick 
gas and ~ 1% is covered by predominantly neutral gas. Note that 
despite the lower covering factor, this gas may still contribute to 
the integrated opacity of the ionizing radiation (see below). About 
30% of the area is occupied by gas with A'^hi ^ 10^^ cm~^. Most 
of the MFP gas and of the LLSs are associated with the streams 
alone, and half of the area that is covered by SLLSs resides outside 
the virial halo of the massive central galaxy, with streams account- 
ing for two third of the covering factor. A significant fraction of 
the neutral gas is located within halos with half of the DLA cross 
section inside the streams. Above A'^ei ~ 10^' cm~^, the fraction 
of DLAs in the streams drops below 20% and becomes negligible 
at even higher column densities. Only 20% of the DLA cross sec- 
tion is found within 0.2iivir, where it is likely associated with the 
central disk, indicating that a significant fraction of DLAs can be 
found within clumps and satellites inside the virial radius of the 
main halo, but outside the central disk. 

That a non-negligible fraction of the DLA cross section is 



found within the virial radius but outside the central disks appears 
to be a typical property of these simulations, as shown in Figure|5] 
In each snapshot (STAR model), we identify an independent struc- 
ture as an ensemble of contiguous cells above Nui = 2 x lO^'' 
cm~^. We define the central galaxy to be the closest structure to 
the center of the main halo and we consider satellites or clumps in 
the streams to be all of the remaining groups of contiguous cells. 
The average cross section as a function of the halo mass is shown 
with blue squares for the central galaxies and black circles for the 
central galaxies and satellites. The amplitude of the halo mass inter- 
vals, chosen to have an equal number of objects, is indicated by the 
horizontal error bar, while the vertical error bar is for the standard 
deviation along the mean. Values are listed in Table lB6l 

In the entire sample, the central galaxy contributes ~ 35% — 
65% of the DLA cross section, with galaxies residing in more mas- 
sive halos showing smaller fractions. We conclude that in these 
simulations the satellites of massive halos and clumps that re- 
side in the streams are at least as important to the DLA cross- 
section of massive dark matter halos as the central galaxy (see also 
iMaller et al1l200lt iRazoumov et alj|2006h . Note that the number 
of objects included in this study is limited and the cross-sections 
listed here can be affected by sample variance. Also, the cross 
section of DLAs is particularly sensitive to the disk sizes that 
are notoriously difficult to correctly reproduce in simulations, al- 
though our simul ations show extended disks with 3 — 5 kpc radii 
at high-redshift. dCeverino et"Zll2010l) . Further, these simulations 
have been selected to study cold gas accretion in massive galax- 
ies and are not representative of the entire population of galaxies 
giving rise to ALSs. Nevertheless, it is reassuring that the mean 
DLA cross section within _Rvir ranges between 50 — 800 kpc^ for 
halos of 10^" — 10^' Mq, in good agreement with the cross sec- 
tio ns found at similar m asses in the larger cosmol ogical simulation 
bv |Pontzenetal.l ^20081) and lTescari et alj ( |2009|) . Our predictions 
are towards the lower end of the size distribution from ted (l2010h . 
who finds a larger area as a consequence of strong feedback. This 
difference emphasis how various feedback prescriptions affect the 
final cross sections. 

The total area is nearly proportional to the halo mass below 
4 X lO''^ Mq (slope 7 = 1.06 and normalization /3 = 3.16 x 
10^" for a log-log linear regression; dashed line in Figure |5] and 
roughly constant at higher masses. Across the whole available mass 
range, an approximate fit to the points in Figure |5]is given by the 
mass dependence of the area within the virial radius at a fixed time: 
Avir oc JJ^ij. oc M^if . This is consistent with DLAs originating 
from regions where the surface density is above a fixed threshold, 
provided that the distribution of densities within the virial radius 
(or 2i?vir) is self-similar among the halos. 

Note that the relation found in Figure|5]is based on only seven 
independent halos at all times and redshift is partially degenerate 
with the halo mass, since the most massive halos are found only at 
late times. It is therefore useful to inspect values for single galaxies, 
superimposed in this figure at z ~ 2.8 and z ~ 1.9 in three projec- 
tions. While a mild evolution with redshift is visible (see also next 
section), the mean cross sections previously derived appear to be a 
fair representation of this sample as a whole, more so considering 
the large intrinsic scatter from galaxy to galaxy. 



4.2 Redshift evolution of the covering factor 

A qualitative description of the covering factor evolution within 
2i?vir is provided in Figure |6] where we show a gallery of pro- 
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Figure 6. Gallery of projected H I column density (STAR model) in selected galaxies and redshifts. Four intervals of column density are marked with different 
colors (DLAs in red, SLLSs in green, LLSs in blue and MFP gas in orange, similar to Figure|4)- Redshift is decreasing from top to bottom, and viiial mass is 
increasing at constant redshift from left to right. The dotted circles mark the virial radius. Cold streams are "patchy", with neutral pockets of gas embedded 
in a more widespread ionized medium. A difference is seen between 2 > 3 and z < 2, and a marginal mass dependence is hinted, but the latter could be an 
artifact of numerical resolution. 
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Figure 7. Time evolution of the cumulative covering factor within 2/?vir- DLAs are in red, SLLSs in green, LLSs in blue, and MFP gas in orange. The central 
colunms are for the galaxies and streams in the STAR model, the left columns are for the streams alone in the STAR model and the right columns are for the 
UVB model, galaxies and streams. The numbers above the bars indicate the number of galaxies in each bin of redshift. Gray dashed and solid lines indicate 
the standard deviation for the mean covering factor of DLAs and LLSs respectively. The covering factor slowly decreases with redshift at all column densities, 
reflecting the decline of cosmological density. 



jected H I column densities (STAR model) for a selected number of 
galaxies. Redshift decreases from the top to the bottom and mass at 
fixed redshift increases from the left to the right. DLAs are marked 
in red, SLLSs in green, LLSs in blue and MFP gas is in orange. 
In each panel, the galaxy center and virial radius are indicated by 
a dotted circle. Gas that surrounds the galaxies is non-uniform and 
the cold streams are "patchy" with neutral regions that are embed- 
ded in a widespread ionized medium. 

An evolution in redshift appears from this figure. The gas cov- 
erage is larger at early times than at late times. Evolution in the 
gas concentration is evident for DLAs. At highest redshifts, there 
are large regions of neutral gas within the streams. These become 
more confined to the central disks and satellites at later times. A 
trend in mass is also visible, although less pronounced, with larger 
gas coverage in lower mass halos. Also, the topology of the streams 
seems to evolve with time and mass. In low mass systems and at 
early times, wide streams penetrate inside the virial radius to de- 
liver cold gas to the forming galaxies. At lower redshifts and in 
more massive systems, the streams break into a series of narrow 
filaments (see[K eres & Hemquist 2009; Faucher-Giguere & Keres 
l201lHstewart et al.l20100 . In our simulations, numerical resolution 
could be responsible for the increased smoothness in the smallest 
galaxies at the highest redshifts and higher resolution is required to 
confirm whether this trend is real. 

A more quantitative analysis of the redshift dependence of the 
covering factor is given in Figure |7] Each bar, represents the mean 
values for DLAs (red), SLLSs (green), LLSs (blue) and MFP gas 



(orange) at one redshift and within 27?vir- The central column is for 
the galaxies and streams in the STAR model, the left column is for 
the streams alone in the STAR model and the right column is for the 
UVB model (streams and galaxies). Values for DLAs and LLSs in 
individual galaxies are listed in Table lB4l All covering factors are 
given as cumulative, i.e. integrated from the lowest column density 
to infinity (numerically approximated to A'^hi ~ 10^^ cm~^). The 
numbers at the top of the bars indicate the number of galaxies in- 
cluded in each redshift bin. The standard deviation of the sample is 
shown for DLAs and LLSs with gray dashed or solid lines to high- 
light the wide scatter from galaxy to galaxy, partially associated to 
the mass (see Figure |9). 

Figure [7] confirms the impression that the covering factor 
within 2_Rvir is slowly declining with time. At the highest redshifts, 
30 — 50% of the area is covered by MFP gas, 7 — 10% is cov- 
ered by LLSs and 1 — 2% is occupied by DLAs (compare with 
iFaucher-Giguere & Keresll20r ij). By z ~ 2, the covering factor be- 
comes 1% for DLAs, 5% for LLSs and ~ 20% for MFP gas. In 
the range z = 2 — 1.3, the covering factors are ~ 10%, 3 — 4%, 
and 0.3 — 0.4% respectively. Note that the already limited statistics 
are dominated by lower mass galaxies below z ~ 1.9, while the 
highest redshift bins may be more prone to resolution effects. 

In the redshift interval 2 < z < 3, where the sample is more 
homogeneous, the time evolution of the covering factor can be 
parametrized by ~ o"'', with a the expansion factor. This trend can 
be explained by the universal expansion for a self-similar density 
distribution of neutral hydrogen. In these simulations, the cumu- 
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Figure 8. Same as Figure [T] comparing covering factors that are measured in different areas. In the center, the covering factor is measured within a fixed 
aperture of 10 arcsec in radius. Shown for comparisons are the covering factor defined within /?vir and 2i?vir> on the right and left, respectively. The number 
of galaxies included in each bin is shown above the columns (variable/fixed apertures) and the standard deviation around the mean is highlighted with solid and 
dashed lines. The variation with time is much weaker when the aperture is fixed, reflecting the decrease of covering area with distance from the halo center 
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Figure 9. Redshift evolution of the covering factor for LLSs in individual 
galaxies (color coded as in Figure[T) within 2iJvir (bottom panel) and for a 
fixed aperture of 10 arcsec (top panel). Trends similar to those infeiTed from 
Figure |7] and [8] are visible when inspecting individual galaxies, regardless 
to their mass. 



lative Hi column density PDF can be described by oc A^j^i'^, with 
c ~ 0.2 for A'hi ^ 10^° cm~^ and c ~ 1 at higher column den- 
sities. In our sample, the viral mass and the virial radius grow with 
time as ~ a? and ~ a?, and the virial column density decreases as 



~ a ^ . For a self-similar density distribution of neutral hydrogen, 
it follows that c ~ d, in agreement with the evolution in Figure|7] 

Comparing next the UVB model (right bars) with the STAR 
model, a systematic offset towards higher values is apparent from 
the former. The inclusion of local sources does not drastically al- 
ter the covering factor (a factor of two for DLAs and LLSs), but 
it is clear that ionizing radiation from stars has appreciable effects 
especially at intermediate column densiti es and it is req uired for a 
consistent comparison with observations JSchavel [20061) . The red- 
shift evolution is still present in the UVB model, reassuring us that 
it is not driven by local sources. Furthermore, this trend cannot be 
artificially induced by our choice of a constant UVB, since the in- 
tensity of the UVB decreases towards the highest redshifts here 
considered. 

In Figure [7] we finally compare the covering factor from the 
streams and galaxies with the streams alone. Overall, the cover- 
ing factor of highly optically thick gas in the streams is decreasing 
faster with time than that of galaxies. For the MFP gas, this effect is 
absent, with streams responsible for the totality of the covering fac- 
tor at all times. Streams contribute to more than 80% of the LLSs 
covering factor over the full redshift interval. Instead, a rapid evo- 
lution is seen for SLLSs and DLAs. Streams are responsible for 
as much as 70% of the covering factor at the highest redshifts, but 
less then 30% of the covering factor at later times as neutral gas 
accumulates in proximity of the central galaxies an satellites. 

Being a normalized quantity, the covering factor within 2i?vir 
is suitable for comparisons among different galaxies and redshifts. 
However, it does not provide direct information on the physical 
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cross section and its evolution. Further, observations probe either 
the gas cross section or the covering factor within a given projected 
distance from a source. In the latter case, the most natural quantity 
to compute is the covering factor within a fixed angular distance 
from the halo center. In Figure |8] we show the mean covering fac- 
tor within fJvir (right columns), 2_Rvir (left columns) and within 
a circular aperture of 10 arcsec in radius (central column). Again, 
the number at the top indicates the number of galaxies included in 
each bin with the second number for the fixed aperture. The vari- 
ance about the mean for LLSs and DLAs is shown with a solid 
and dashed line, respectively. At the redshifts examined here, an 
angular distance of 10 arcsec coiTesponds to 70 — 85 proper kpc, 
equivalent to ~ 3 — 5 times the virial radius at the highest redshifts 
and to ~ 0.5 — 1 times the virial radius at the lowest redshifts. The 
snapshots in which the extraction box exceeds the simulated region 
are not included. Further, for a consistent comparison, we consider 
a path length through the grid of twice the virial radius in all three 
cases. 

The covering factor measured within _Rvir is higher than the 
one inside 2J?vir, particularly for DLAs. This is trivially expected 
for a gas distribution with declining column density. Considering 
now a fixed angular distance, the covering factor exhibits a much 
weaker evolution with a modest increase in the interval 2 < 2 < 3, 
where the sample is more homogeneous. As the universe expands, 
the mean column density within 2i?vir drops and the fraction of 
the area within 2_Rvir subtended by gas above a fixed column den- 
sity threshold decreases accordingly. However, the virial radius is 
increasing with respect to the size of the extraction aperture, result- 
ing in a nearly constant or slightly increasing covering factor. 

Within ~ 100/ft kpc, the covering factor of absorbers in 
the transverse direction of QSO host galaxies is significantly 
higher than our estimates, b eing unity for A'^hi > 1-6 x 10^^ 



cm taennawi et alj l2006h and 30% for A'^hi > lO" cm"^ 



jHennawi & Prochaska l20()7h . This discrepancy may originate in 
the environment around QSOs, more massive than the galaxies in 
our sample. Further, these simulations do not resolve overdensi- 
ties on very small scales and consequently the cross section dis- 
tribution suffers from an additional uncertainty. For example, a 
population of neutral clouds of < 100 pc in size that are p res- 
sure confined by the hot halo gas (see e.g. D ekel & Birn boiml 
I2OO8I : IProchaska & Hennawill2009l ; iBimboim & Dekelll201 ih may 
be missed in these simulations and the cross section of optically 
thick gas would be underestimated. Similarly, shells or bubbles as- 
sociated with radiative shocks on small scales may not be properly 
resolved (cfr. Simcoe et al. 2006). 

The weak redshift evolution seen in Figure |7] and Figure |8] 
may be the result of having combined galaxies at different mass. 
To explicitly test this hypothesis, in Figure |9] we display the LLS 
covering factor within 2i?vir (bottom panel) and within 10 arcsec 
(top panel) for individual galaxies. Although with significant scat- 
ter from object to object, the covering factor within 2i?vir appears 
to be roughly constant or slowly declining regardless to the galaxy 
mass. Similarly, the covering factor within a 10 arcsec aperture is 
slowly increasing at all masses. 



factors previously measured into more direct observables that can 
be used in comparison to real data to uncover signatures of cold gas 
accretion. We start by introducing some of the formalism adopted 
in typical measurements of ALSs. Readers familiar with this nota- 
tion can skip the next sub-section and continue from Section [J!2l 
where we derive the column density distribution function for cold- 
stream fed galaxies and discuss their visibility in absorption. 

5.1 Common formalism 

To probe any form of gas that lies along the line of sight, observers 
collect large spectroscopic samples of QS Os that serve as ran- 
domly distributed ba ckground sources (e.g. IProchaska et alj|2005l : 
lO'Meara et alj2007h . In this experiment, ther e is no a priori knowl- 
edge of the foreground system (see however iRafelski et al.ll2O10l) 
and one can directly measure the probability to intersect gas along 
a random sightline. This quantity translates to a sky covering factor 
that is the percentage of total observed area covered by gas. 
One can define the incidence 

l{z) = ^ (3) 

as the number of systems A^abs detected across the total searched 
redshift Az. This quantity is related to the probability of intersect- 
ing a system along the line of sight at redshift z. In analogy with 
the idealized experiment of a body moving in a sea of particles with 
number density n and cross section cr, the following equality holds 



£{z)dz = nabsCTabsd'c 



(4) 



with Tiabs the comoving number density of absorbers, CTabs their 
comoving cross section and die the infinitesimal comoving path 
length. 

Within a given cosmology H (2:) = Ho[flA + {1 + z)^Qm]^/^ 
and die = c/H{z) dz. Replacing z by an "absorption length" X 



dX = ^^{1 + zfdz 
M [z) 



(5) 



and using the identity £{z)dz = £{X)dX, the incidence can 
be rewritten as the num ber of systems per absorption length 
( lBahcall&Peebleslll969h 



£(X) = -r--nabs'^abs 
-Ho 



(6) 



where the cross section <^abs now represents the physical area of 
the absorber. 

Equation ^ is the zeroth moment of the column density 
distribution function f{Niii,X) that characterizes the number of 
ALSs per absorption length as a function of the H I column density. 
From the definition of £{X), it follows that 



f{Nm,X) = —n,i,,(Nm,X)A,UNm,X) 
Ho 



(7) 



with j4abs the physical area per unit A'^h i- Equation ^ reveals that 
absorption line statistics directly probe the shape and the evolution 
of the physical cross section times the number density of absorbers 
rather than the covering factor of gas around these systems. 



5 PREDICTED ABSORPTION LINE STATISTICS 

In the previous section we discussed the redshift evolution of the 
covering factor in streams and galaxies, where we found that only 
a small fraction of the virial area is covered by optically thick gas. 
In this section, we wish to translate the cross sections and covering 



5.2 The column density distribution function 

To understand the nature of the gas probed in absorption, it is com- 
mon practice to reproduce the f{Nui, X) in its shape and normal- 
ization from simulations of large volumes, required for statistical 
significance. Here we take a different approach where we derive 
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Figure 10. Comparison of the shape of the simulated and observed column 
density function /{Nn i,X). The normahzations of the different curve are 
arbitrary, to highlight the shapes only. The simulated /{Nm, X) from the 
STAR model at z ~ 3.5 is shown for gas inside the virial radius (solid 
line), for gas within 2_Rvir (dash-dotted line), and with for the streams only 
(dashe d Hne). The observed f{Niii, X) at 2 ~ 3.7 from lProchaska et al.l 
)2O10i) is shown as a colored curve. Power law regressions to the simulation 
results in intervals of column density are shown as dotted lines. The num- 
bers denote the power law indexes. Simulations of massive galaxies at high 
resolution reproduce the observed cutoff at high column density, but predict 
an f{Niii, X) that is shallower than observed around A^hi ~ 10^" cm~^ 
andforSLLSs. Cold streams dominate the /(TVhi,^) below A^Hi ~ 10^* 
cm~^, but a larger contribution from the intergalactic medium is required 
at lower A^h i to satisfy the constraints imposed by the observed mean free 
path of ionizing radiation. 



an f(Nni, X) that, by constmction, is due to gas associated only 
with massive galaxies and the streams with merging galaxies that 
feed them. Any similarity or discrepancy with the observed distri- 
bution function informs the features of f{Nui, X) that reflect the 
gas distribution in these particular class of systems as predicted by 
these simulations. The reader should keep in mind that our sample 
and our f{Nui, X) are not intended to be representative of the full 
ALS population. 

Our model for f{Niii,X) follows the definition in Equation 
Q. The physical area per unit A'h i is constructed from the average 
cross section distribution function over all the galaxies at a common 
redshift 



Agal(iVHI,2;) 



0gal.i(A^HI, z) 
iVgalAiVHI 



(8) 



We consider the virial radius to define the gas that is physically as- 
sociated with a galaxy, but we construct also a second model con- 
sidering the cross section within 2_Rvir to stress the effects of the 
cold streams. For the comoving number density of absorbers, we 
use the cumulative number densit}|3 of galaxies that are more mas- 
sive than the least massive galaxy at a given redshift. In practice, we 
construct an average cross section that we attribute to all the halos 
of equal mass or more massive than the galaxies in our sample. This 



The como ving number density is com puted following the approximations 
described in Dekel & BimboimI ilQOd Appendix A). Recent simulations 
(e.g. Klvpin et al. 201(f reveal deviations from these approximations, but 
they become relevant only above z ~ 4. 



procedure implies an extrapolation towards higher masses, since 
systems much above 10^^ M© are not represented in our sample. 
Due to the scatter from galaxy to galaxy in the cross section and 
to its mass dependence, Ags_i{Nm, z) and the final f{Nui,X) 
should be regarded as approximations. 

In Figure [To] we compare the functional form of the 
f{Niii,X) derived for the STAR model at z ~ 3.5 within 
the virial radius (solid black line), within twice the virial ra- 
dius (dashed-dotted line), and in the streams alo ne (dashed line). 
The o bserved distribution al z ~ 3.7 from IProchaska et al.l 
( I2OIOI) is superimposed with colors. For DLAs and SLLSs (i.e. 
A'h! > 10^' cm~^ where one resolves the damping wings of 
Lyof), f{Niii,X) is derived from observations of individual ab- 
sorbers, but at lower column densities observations provide only 
integrated constraints on the distribution function. For this rea- 
son, the shape of f{NHi,X) in the LLSs and MFP region is 
more uncertain. In Fig ure 110] we report the preferred model from 
IProchaska etal.l ( l201(j) . 

Although not a perfect match, these models reproduce some of 
the features of the observed f{Niii,X). Three different intervals 
can be identified. First, there is a steep decrease at high column 
density with an evident break around A^hi = 3 x 10^^ cm^^. From 
our model, we measure a power law ~ — 2.7 for log A^hi > 21.7 
(in cm~^; dotted line), similar to the observed slope ~ —3. This 
cutoff naturally emerges from our high resolution simulations of 
massive disks, while it is hard t o reproduce in cosmological sim- 
ulations at lower resoluti on (e.g. Pontzen et al.ll2008: Tescari et al] 
l2OO9l : ICeiill2O10l, but see lAItav et al]|20ld ; lMcOuinn et al.ll201ll) . 
The inclusion of molecules would suppress even further the num- 
ber of systems with Hi above 10'^^ cm~^(see Figure IaTT i. bringing 
the model to even better agreement with observations. 

Moving towards lower column densities, the model exhibits 
a plateau with slope ~ -1 between A^hi = 3 x 10^* - 10^^ 
cm~^. The observed f{Niii,X) distribution also flattens below 
lO^^cm"^ (O'Meara et al. 2007; Prochaska et al. 2010) which has 
been interpreted as the result of transiti oning from predominantly 
neutr a l gas to predominan tly ionized dZheng & Miralda-Escudel 
I2OO2I : iMcOuinn etalj|201lh . In fact, our model flattens at some- 
what high A'h! and underpredicts the observed incidence of SLLSs 
relative to DLAs. Provided that SLLSs arises from galaxies, per- 
haps this discrepancy would be resolved by including less massive 
galaxies than those modeled here. As previously noted for the cross 
section, we might be missing a population of small clouds of neu- 
tral hydrogen due to limited resolution. Finally, below Nni = 10^* 
cm^^ the distribution function becomes moderately steeper with 
slope ~ —1.3. This qualitatively follows the observed distribution, 
but it is not sufficiently steep to reprod uce the integrated mean free 
path to ionizing radiation at 2: ~ 4 jProchaska et al.ll2009l) . We 
infer, therefore, that gas beyond the virial radiu s must contribute 
signifi cantly at these column densities (see also lKohler & GnedinI 
l2007h . As already seen in Figure |4] due to the stream contribution 
(dashed line) the f(NHi,X) steepens mildly below A^'hi ~ 10^* 
cm^^ , but gas within 2_Rvir is not enough to reproduce the observed 
slope (dash dotted line). Contribution from the IGM at larger radii 
may be required. 

Figure [TTI compares again the observed /(7Vhi,X) with the 
model inside the virial radius, including at this time the abso- 
lute normalization. The observed f{Nui,X) has been corrected 
to match our adopted cosmological parameters. Values are listed 
in Table |B5] This figure reveals an overall satisfactory agreement 
between the observations and the models for high column den- 
sity DLAs and at A'^hi ~ 10^* cm~^. As already noted, the two 
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Figure 11. Column density distribution function from observations at 2 '~ 3.7 f colors : |Prochaska et alJllOlol) and simulations for the STAR (solid line) and 
the UVB (dashed line) models at 2 = 3.5 within the virial radius. Shown in gray are simulated /(A^jj i, X) for different redshifts in the interval 2 = 4 — 1.4. 
There is a satisfactory agreement between theory and observations for DLAs and at A'jji ~ 10^* cm"'^, while the two distributions diverge at lower column 
densities. SLLSs and high density LLSs are underrepresented by these models. The /(A'jji, X) exhibits modest redshift evolution in shape, in agreement 
with observations. 



distributions diverge entering the optically thin regime, suggesting 
a missing IGM contribution in the simulations. Also missing are 
SLLSs and high column density LLSs that are not found in large 
enough numbers within the halos of massive galaxies. The UVB 
model, superimposed as a dashed line, has a larger fraction of LLSs 
and SLLSs, but not enough to account for the missing systems. It 
is interesting to note that in these models SLLSs originate from gas 
at the outskirts of the disks and are ionized b y the leaking U V pho- 
tons from stars that are forming in the disk ('Schave"2006 |). While 
most of the f(NB_i,X) is shaped by the UVB ( Razoun iov et al.l 
l2006l : lNagamlne et alJlOlO ). RT effects from local sources produce 
non-negligible effects. Albeit with all the uncertainties discussed, 
from this simple comparison we infer that massive disks shape the 
/(A^Hi, X) at the highest column densities and cold streams con- 
tribute non negligibly to the LLS population, particularly between 
A^Hi ~ lO^'^ - 10^** cm^^ 

Finally, in gray, we superimpose different f{NHi,X) for 
the redshift interval z = 4—1.4. These models are remark- 
ably invariant with redshift and this limited evolution in the 
/{Nh is consistent with the observation that the column den- 
sity distribution function o f DLAs preserves its shape over 12 Gyr 
IProchaska & Wolfell2009l) . This is also consistent with the present 
day distribution function from 21-cm observations in H I rich disks. 
Approaching t he optically thin reg i me, a larger degree of evolution 
is observed (iProchaska et al ]|2010l.lRibaudo et al.l201 lal) . which is 
attributed to radiative transfer effects. This variation is not captured 
by our model, probably because massive galaxies are less sensitive 



to the effect of the UVB than the IGM and low mass galaxies which 
are needed in order to fully match the observed distribution. 

The limited evolution in the models presented in Figure [TTIcan 
be qualitatively understood from visual inspection of Figure |6] As 
the virial radius grows with time, the gas preserves its distribution 
with centrally concentrated neutral regions surrounded by SLLSs. 
In turn, this column density interval marks the transition to more 
filamentary structures that dominate the LLSs and MFP cross sec- 
tion. More quantitatively, this non-evolution is related to the fact 
that the neutral hydrogen column density PDF is found to be nearly 
self-similar. Therefore, in our simulations, the lack of evolution in 
f{NHi,X) is driven by the distribution of gas within each halo (i.e. 
the cross section) more than by the cosmological evolution (i.e. the 
number density of systems). 



5.3 Are cold streams visible in absorption? 

By integrating the column density distribution function, one can 
compare predictions from the simulations with the observed inci- 
dence of absorbers. Figure [TT] shows that the normalization of the 
simulated /(Nui, X) is not too far from the observed one. In these 
simulations, massive galaxies (Mvir ~ 10^'' - 10^^ Mq) alone can 
account for ~ 20% — 30% of the observed DLA and LLS popu- 
lation. The inclusion of satellites and streams up to twice the viria l 
radius further increases this contribution (cfr. iMaller et alj|2003l) . 
Although uncertain up to a factor two, these numbers suggest that a 
non-negligible fraction of observed ALSs originate in galaxies and 
streams as represented in these simulations. 



Absorption line systems in simulated galaxies fed by cold streams 15 



This result is not at odds witli tiie low covering factor for op- 
tically thick gas discussed in the previous sections. This is because 
the incidence of absorbers with background QSOs is sensitive to 
the physical cross section times the number density of halos. Given 
a large enough number density of halos, even galaxies with rela- 
tively small covering factor can offer enough total cross section to 
account for a significant fraction of the observed incidence. Most 
of the DLAs are estimate d by theory to originate in the mass inter- 
val 10^ - 10" Mq (e.g. lPontzen et alllloOS ; Tescari et al. 2009i). 
with a pe ak around ^ 5 x 10^° Mp. In m odels with strong feed- 
back (e.g. iNagamine et alj200'7l : lcej2010l) . the peak shifts to even 
higher masses. It is not surprising that our sample that includes 
galaxies with 5 x 10^" — 3 x 10^^ M© at 2: ~ 3 is within a factor 
of a few of the observed incidence. 

Combining the observed incidence for DLAs and LLSs with 
the cross section from our simulations, we can infer the mini- 
mum halo mass required to fully match observations. At 2: ~ 3.5 
{z ~ 2.8), the average cross section within 2i?vir is ~ 304 kpc^ 
(~ 530 kpc^) for DLAs and ~ 1548 kpc^ (~ 2010 kpc^) 
for LLSs. The halo number density r equired to match the ob- 
served incidence of DLAs and LLSs (Procha ska & Wolfell2009l; 
[Prochaska et al. 2010; Ribaudo et al 2011a) is ~ 0.062 Mpc"^ 
(~ 0.032 Mpc^''), corresponding to a minimum halo mass of 
1.6 X 10^° M0 (3.6 X 10^'' M0). This is a factor 1.4 (2.7) be- 
low the minimum mass included in this sample. This calculation is 
very crude since, for example, at lower masses the mean cross sec- 
tion is expected to decrease (e.g. Figure |5}- Further, these results 
are sensitive to the radius adopted for the cross section determina- 
tion, here chosen to be 2_Rvir, i.e. the size of the re-simulated box. 
Since ideally one would adopt the mean separation between halos, 
larger than 2_Rvir at these masses, our numbers may be even un- 
derestimated. With these limitations, our calculation suggests that 
only a small extrapolation towards lower masses is needed to fully 
account for the observed incidence of ALSs. 

In summary, despite their small covering factor in each halo, 
the cold streams should have already been detected in absorption 
in large spectroscopic surveys. Since most of the cross section in 
cold streams is predicted to be at A'^ei ;$ 10^'"* cm~^(see Figure 
|4]and Figure [TOt, LLSs are the best candidates for the cold gas in 
the inflowing streams. If there are other components of cold gas 
in the halo that may not exist in proper abundances in our current 
simulations (e.g. high velocity cloud analogs or massive outflows), 
then other indicators are required to distinguish between them and 
the cold streams. 



6 OBSERVABLE PROPERTIES OF THE CGM 

In a second class of experiments designed to probe cold gas in ab- 
soiption, observers select an object or a class of objects of interest 
and use background sources to probe the gas distribution around 
the foreground systems. Using a background quasar (or galaxy) 
and knowing the angular separation between the probe and the ob- 
ject of interest, this experim ent yields a direct measurement of the 
angular Hi distribution (e.g.lChristensen et al1l2007t iMonier et al.l 
I2OO9I : iFumagalli et a lj|2010bf) about the object of interest or of the 
gas covering factor and kinematics (e.g. Steidel et al. 2010). In the 
next section, we provide a summary of the formalism adopted in 
the analysis of absorption lines as well as a more technical descrip- 
tion on how line profiles are derived from these simulations. Start- 
ing from Section|6j2l we discuss predictions for the metallicity and 



kinematics of absorption lines, providing a comparison with recent 
observations. 



6.1 Common formalism and numerical procedures 

6.1.1 The equivalent width 

The normalized profile of an absorption line is characterized by the 
equivalent width 



c 



I{v)] dv , 



(9) 



with I the normalized intensity and Aq the rest frame wavelength at 
the line center. The equivalent width provides a direct measurement 
of the light absorbed along the line of sight. Through the radiative 
transfer equation, W{\()) can be related to the gas optical depth 
that in turn depends on the gas column density and kinematics and, 
via the oscillator strength, on the quantum mechanical properties of 
the transition. 

Three different regimes are typically identified in the so-called 
curve of growth, that describes the variation of the equivalent width 
with both column density and oscillator strength. At very high col- 
umn densities, in the damped regime, all the light at the line cen- 
ter is absorbed and the equivalent width is a solely function of the 
column density. At intermediate column densities, on the saturated 
portion of the curve of growth, the line center is again fully black, 
but the equivalent width is most sensitive to the line of sight com- 
ponent of the gas velocity that regulates the line width. Finally, at 
low column densities, in the linear regime, the equivalent width be- 
comes proportional to the column density in the observed direction 
and the absorption profile is no longer zero at the line-center. Av- 
eraging along different sightlines, one can further probe the mean 
covering factor around a class of sources by looking at the fraction 
of transmitted light close to the line center. 



6.1.2 Techniques for the simulated line profiles 

To reproduce the absorption line-profile for Lya and SillA1260, we 
compute the transmitted intensity Iv of a background source with 
intensity 7o in a random projected sightline i 



= lo exp(— r^,i) , 
where the optical depth r,y,i is given by 

k\i 



(10) 



(11) 



Here, the sum is over all the AMR cells aligned along the i-th line 
of sight, and is the column density in each cell. The frequency 
integrated absorption cross section s is given by 

2 

(12) 



s = / , 

nieC 



with / the oscillator strength, e the electron charge, me the electron 
mass and c the speed of light. The frequency-dependent line-profile 
is 



Hk{u, a) 



with the Voigt function 

H{u,a) = - / dy- 2 

J-00 {u~yy + a^ 



(13) 



(14) 
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Figure 13. Distribution of metallicity, column density weighted, in the en- 
tire sample for DLAs (red dashed hne), SLLSs (green dotted line) and LLSs 
(blue dot-dashed line). Also shown are the streams alone (black solid line) 
and satellites without centrals (orange long dashed line). The smooth stream 
component is more metal poor than the galaxies, with typical Z ~ O.OIZq . 
However, gas in streams is not piistine, but already emiched by previous 
episodes of star formation. 



density is necessary to significantly alter our results. In this simple 
model, we also neglect Si ll recombination at low temperature and 
a lower collisional rate in shielded regions. These effects alter the 
Sill column density, but only in regions of small covering factor 
with minor consequences on the average absorption lines. A com- 
plete analysis of metal lines requires the inclusion of collisional 
ionization and photoionization, which is beyond the scope of this 
paper. 



Figure 12. Metallicity in the optically thick gas of MW3 at 2 = 2.3, 
weighted by hydrogen column density. Top: galaxies and streams. Bottom: 
streams alone. The black and white contours denote log(NHi/cm~^) = 
17.2 and log(NHi/cm~^) = 20.3 respectively. A metallicity gradient is 
visible. Gas in and around galaxies is more enriched than in the streams 
alone. 

Here, a — ^/(AtvAvd), u = {y — vi)/A.vd, y ~ v/b and 
Aud = {uob)/c. In the above equations, 6^ — {2kT/mx + C^) 
is the broadening parameter given by that gas temperature and el- 
ement mass, 7 is the sum over the spontaneous emission coeffi- 
cients, 1/ is the observed spectral frequency, and = 10 km/s the 
turbulent velocity. We assume the central galaxy is at zero velocity 
(coiTesponding to the rest-frame frequency i/q) and we compute for 
each gas cell the line center frequency i/i corrected for the gas ve- 
locity along the line o f sight. We integ rate Equation ( I14t using the 
formalism provided in'Zaghloul' (zOOT). 

All the spectra are computed at a resolution of 4 km/s. For 
Lya, we assume / = 0.4164 and 7 = 6.265 x 10*. For Si nAI260, 
we use / = 1.007, and 7 = 2.533 x 10^ and we derive the Si 11 
density from the neutral hydrog en density and meta llicity, assum- 
ing log(Si/H)0 + 12 = 7.51 ('Asplund et a"i]|2009j) . To compen- 
sate for the different ionization potential of silicon and hydrogen, 
we boost the Si 11 volume density by a factor of 10 in cells with 
xhi < 0.1. The amplitude of this boost depends on the gas and ra- 
diation field properties, but it will become clear from our analysis 
that a variation of several order of magnitudes in the Si 11 column 



6.2 Metallicity 

The metal enrichment of the gas may help distinguish between 
gas that is accreting from the cosmic web and the more enriched 
material that is located in proximity to a galaxy or is outflowing 
from regions with intense star formation. Production and diffu- 
sion of metals is known t o be problematic in numerical simula- 
tions IWadslev et al.ll20o3) . but the AMR code used for our sim- 
ulations should minimize problems of me tal transport typical of 
SPH simulations (but see IShen et alj[201ol) . Caution is advised in 
generalizing these comparisons of metallicity between our simu- 
lations and observations because our current simulations are lim- 
ited to supernova-driven winds which do not reprodu ce outflows 
as ma ssive as suggested by some observations (e.g. Steidel et alj 
|20IO|) . and also because the specific IMF , Miller & Scalo. (1979.) is 
assumed in the simulations. 

In Figure [T2I the column density weighted metallicity in solar 
units (Zq = 0.0134; Asplund et al. 2009) for a projection of MW3 
at z ~ 2.3 is shown. On the top, we display metals in both streams 
and galaxies for optically thick gas; on the bottom, we present the 
metallicity in the streams alone. Contours are for A^hi = 1-6 x 
lO^'^ cm-2(black hne) and Nm = 2 x 10^° cm-2(white hne). A 
spatial gradient in metallicity is visible. DLA gas inside the main 
galaxy and satellites is significantly enriched, with values close to 
1/10 solar (in r easonable agreement with DLA observations; e.g. 
IProchaska et al.l 12003). The fraction of metals decreases moving 
far from the central disks and at lower column densities. The bulk 
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of the LLSs is enriched at 1/10 — 1/100 solar, with the streams 
skewed at the lowest metallicity. 

A similar metal content characterizes the entire sample at 
z = 1.3 — 4, almost independently of redshift and mass. Figure 
[T3] shows the distribution of column density weighted metallicity 
tp in DLAs (red dashed line), SLLSs (green dotted lines) and LLSs 
(blue dash dotted line). The metal distributions for streams alone 
(only for A'^ei > 1-6 x 10^^ cm~^; solid line) and for satellites 
without centrals (orange long-dashed line) are superimposed. This 
distribution has been normalized so that the integral over the metal- 
licity gives the average covering factor in the entire sample for a 
fixed range of column density 



ijiAN-ai,Z) 



<t>{AN-ai, AZ) 
Az (htot 



(15) 



Here, (p is the physical cross section. 

From Figure[T3] we see that most of the area that is covered by 
optically thick gas is occupied by metal poor cold filaments, with a 
distribution centered at Z ~ 10~^Zq and a spread of over ~ 1 dex. 
Cold streams are responsible for a metal poor population of LLSs, 
while optically thick absorbers that originates within 0.257?vir of 
centrals and satellites can be enriched above log Z/Zq = —1. 
DLAs are characterized by a broad distribution of metals centered 
at log Z/Zq ~ — 1 with a tail to lower metallicity. Satellites alone 
have a distribution similar to the one of DLAs. SLLSs originate 
from ionized material in the surroundings of the disks and have 
an equally broad distribution, centered around log Z/Zq ~ —1.5. 
Outflows more vigorous than the ones generated in our simulations 
are expected to produce even higher metallicity around the star- 
forming regions. Since more primordial gas that replenishes galaxy 
disks with fresh fuel for star formation is predicted to have one or- 
der of magnitude less metals than gas in galaxies, we conclude that 
metal poor LLSs in proximity to galaxies are the best candidates 
for cold g as in inflowing strea ms. The metallicity distribution in 
LLSs (e.g. |ProchteretaLll2010l) may offer a crucial test to distin- 
guish between cold inflowing streams and gas that already resides 
in the galaxies or is outflowing (see also .Giavalisco et al..,201 1,) . 



6.3 Absorption line profiles 

The information on the metallicity and covering factor can be com- 
bined together with kinematics to formulate predictions on the 
strength and shape of absorption line profiles at different projected 
distances from the galaxy centers. These models can then be com- 
pared with the population of AL Ss observed in the CGM of Ly- 
man break galaxies at z ~ 2 — 3 jSteidel et al.ll201(il) . We perform 
this calculation by stacking several spectra generated along multi- 
ple sightlines in a subset of our sample (MWl, MW2, MW3 and 
MW8 in two projections) at 2 ~ 2.3 and z ~ 3.2. The absorp- 
tion profiles at z ~ 2.3 for Lya (top panel) and SillA1260 (bottom 
panel) are shown in Figure [T4l where each spectrum is a compos- 
ite of independent sightlines in five intervals of impact parameter, 
chosen to match the observations presented in Steidel et al. (2010). 

Despite the fact that Lya saturates at low Hi column densi- 
ties, none of the spectra are black at the line-center due to non- 
unity covering factor of neutral gas. The transmitted normalized 
intensity at the line center ranges in the interval 0.5 — 0.6, similar 
to what is observed. Damping wings are visible within the inner- 
most 40 kpc, where most of the highly optically thick gas resides. 
The line profiles are characterized by a full width at half maximum 
(FWHM) velocity of ~ 400 - 600 km/s for 6 < 41 kpc, related 
to the kinematics of the central galaxies and the inner satellite rela- 
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Figure 14. Averaged Lya (top) and Si IIA1260 (bottom) absorption lines for 
background sources at different impact parameters from the central galaxy. 
The LyQ profiles are never saturated due to the low covering factor and ai'e 
characterized by a FWHM of ~ 200 - 300 km/s for b > 41 kpc. The low 
covering factor, low intrinsic Si/H abundance and low metallicity produce 
very weak average Si IIA1260 absorption. 



five motions, as well as the bulk velocity of the incoming streams. 
Conversely, at larger impact parameters (6 > 41 kpc), a velocity 
of 200 — 300 km/s FWHM is most likely associated with the in- 
flowin g gas that is streaming at ~ 200 km/s or more dDekel et al.l 
l2009bl) . The line profiles are symmetric, despite the fact that the gas 
is infalling towards the central galaxies, indicating that there are no 
preferred directions in the velocity field in a stack of multiple sight- 
lines. Inspecting the metal lines, it is evident that the low covering 
factor and the intrinsic low Si/H abundance, together with the low 
metallicity of the streams, produce absorption lines that are weaker 
than the observed values at all impact parameters 6 > 10 kpc. 

In Figure \T5\ we present a comparison between the line 
strength from simulations and observations. For each radial bin, 
the rest frame equivalent width of Lya at z ~ 2.3 and z ^ 3.2 
are shown with red open squares and crosses, respectively. Simi- 
larly, we superimpose the equivalent width of Si llA1260 at z ~ 2.3 
(blue open circles) and z ~ 3.2 (blue x's). Values are listed in Ta- 
ble lBVI The black dotted and dashed lines are for the observed rest 
frame equivalent w idth from spectra of LBGs at redshift 2 ~ 2 — 3 
dSteidel et J]|201()h . As expected from direct inspection of the ab- 
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Figure 15. Equivalent width profile al z ^ 2.3 and z ~ 3.2 for Lya 
(red open squai'es and crosses) and SiIlA1260 (blue open circles and x's). 
The black dotted and dashed lines show the observed rest frame equivalent 
widths in LBGs at redshift z ~ 2-3 I Steidel et alj2 010). Massive galax- 
ies that are accreting gas through cold streams have a large enough covering 
factor and line-of-sight velocity to account for the observed strength of the 
Lyo absorption. Conversely, the metal poor streams cannot reproduce the 
observed equivalent width in metal lines. Boosting the degree of enrichment 
by introducing an artificial metallicity floor at O.IZq (green downward tri- 
angles) is not enough for matching the observations. 
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Figure 16. Optical depth integrated 

±250 km/s for SiIlA1260 relative to Lya for gas with A^hi > 10 
cm~^. The low intrinsic Si/H abundance and metal content of cold streams 
suppress the metal line opacity by more than five orders of magnitude com- 
pared to Lya. 



soiption profiles, massive galaxies that are accreting gas througli 
cold streams have a large enough covering factor and sufficiently 
complex kinematics in Hi gas to match observed Lya values. We 
emphasize that this is achieved through the gas associated with 
satellites and streams alone; the supernova-driven winds in our 
current si mulations are hot an d not as strong as the outflows sug- 
gested bv lSteidel et Zl ( I2OI0I) . While winds, commonly observed 
in star-forming galaxies, are essential ingredient for realistic mod- 
els, the fact that the same observations can be reproduced by two 
extreme scenarios stands as a reminder that evidences for outflow- 



ing/inflowing gas are subtle and more refined models are now re- 
quired. 

That these simulations fail to reproduce the strength of Si 11 at 
all impact parameters fe < 100 kpc is not surprising since, due to the 
lack of strong outflows of cold gas, low-ion metals are tracing the 
distribution of hydrogen in the galaxies and streams. We can con- 
sider whether, for example, the discrepancy for SillA1260 might 
be attributable to the low enrichment of the cold streams. In Figure 
[T6l we show the ratio of the Si 11 to Lya optical depth, integrated 
in the velocity interval ±250 km/s. Only gas with A^hi > lO^'' 
cm~^ is displayed. The low abundance of silicon relative to hydro- 
gen at solar metallicity suppresses the metal line opacity by five 
orders of magnitude and the low metallicity of cold streams de- 
creases the Sill opacity even further, by a factor of ~ 10 — 100. 
For these reasons, metal lines that originate in cold streams are 
weak and likely to be undetected. A higher degree of etvrichment 
is not enough to boost the equivalent width by one ord er of magni- 
tude, as demonstrated by a simple experiment (see also lKimm et al.l 
I20T j). If we artificially impose a floor in metallicity at Z = 0.1 Zq 
for all the gas with A'^hi > 10^^ cm~^, we obtain the equivalent 
widths shown by the green downward triangles in Figure [15] A 
much higher degree of enrichment increases the line strength only 
by a modest value. 

We conclude that cold streams are unlikely to produce the 
large equivalent widths of low-ion metal absorption ar ound massive 
galaxies. If the average equivalent widths reported bv lSteidel et al.l 
(2010) are confirmed (e.g. with higher resolution observations), an 
alternate source of opacity must be included in the simulations. 
This could include, for example, absorption by neighboring galax- 
ies (not modeled here), small clumps of cold gas embedded in a 
hotter medium, and/or galactic-scale outflows. However, in order to 
satisfy the constraints imposed by the Lya absorption, these phe- 
nomena cannot carry a very large column density of neutral hy- 
drogen nor contribute to the kinematics apparent in the cold flows 
unless they replace the streams altogether This would fail to supply 
the gas necessary to produce the observed high SFR plus a similar 
outflow rate. Unless these SFRs and outflows are sustained by a 
large amount of reservoir gas, the inflows and outflows must live 
together. This is a natural consequence of the fact that the inflows 
are along narrow dense streams, while the outflows are covering 
a large solid angle between the streams. Additional work in sim- 
ulations is needed to explore stronger wind models and to under- 
stand the interaction between outflowing and inflowing gas (e.g. 
lRazoumovll2009l : |Powell et alj|201(]| ; lvan de Voort et al.ir201ll) . be- 
fore we can reconcile both the Lya and metal absorption profiles 
with a combination of inflows and outflows. 



7 SUMMARY AND CONCLUSIONS 

We investigated the characteristics of cold gas in cold-stream fed 
galaxies in absorption using a sample of seven galaxies drawn at 
random for zoom-in cosmological hydro-AMR simulations with 
high resolution (35 — 70 pc), in the halo mass range 10^'^ — 10^^ 
M0 between 2; ~ 1.4 — 4. We considered the contributions of the 
cold gas in the central galaxies and of the inflowing cold gas in the 
two components of the streams that feed the high-redshift massive 
galaxies from the cosmic web — the incoming small galaxies and 
the larger smooth component. Our analysis focused on the mass de- 
pendence and time evolution of optically thick and thin gas cross 
section and covering factor, on the metallicity distribution in the 
streams and massive galaxies, and on the Lya and metal line kine- 
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matics. The simulations were compared with available observations 
of absorption line systems in the foreground of quasars and with 
kinematics of the circumgalactic medium in Lyman break galaxies 
to highlight what observables are most useful to detect cold streams 
in absorption. The limited sample size and the fact that very strong 
outflows are not included in these simulations make our predictions 
rather specific to the study of cold streams, and caution is advised 
in generalizing our findings to the fuU population of ALSs. Our 
results and findings can be summarized as follows. 

• After post processing the simulations with a radiative transfer 
code that includes dust, coUisional ionization, and photoionization 
from the UV background and local stellar sources, cold streams ap- 
pear highly ionized and not entirely self-shielded structures. Pock- 
ets of high coluirm density neutral hydrogen, mostly associated 
with the central galaxy and instreaming galaxies, are embedded in 
a ionized medium with Nhi < 10^^ cm~^. The densest regions are 
mostly affected by the UV radiation from local sources because 
the UV background is completely shielded in regions of uh > 
0.1 crri^''. While the UV from stars is largely confined by hydro- 
gen and dust absorption in the star-forming regions, the non-zero 
escape fraction is responsible for the ionization of the gas in the 
immediate surroundings of disks (10^*^ < iVm/crri^^ < 10'^"). 
Since most of the mass is found at the highest column densities, 
UV from stars decreases the total neutral hydrogen mass in galax- 
ies. 

• The smooth component of the cold streams dominates the 
cross section of neutral hydrogen below A^hi = 10^*^ cm^^ within 
IRvir- At 2 < 2, less than half of the cross section in the range 
Nm = 10^^ - 10^° cm"^ is due to the streams. At these low red- 
shifts, most of the cross section for A^hi > 10^° cm~^ is in galax- 
ies, with the smooth component of the streams accounting only for 
a small fraction of the total budget. However, a rapid evolution is 
seen with redshift, with the gas clvmips in the streams responsible 
for more than half of the DLA/SLLS covering factor at « ~ 4. 

• No more than 60% of the DLA cross section is associated with 
the central galaxies. Satellites and gas clumps in streams are im- 
portant or even the dominant contributors to the neutral gas cross 
section, especially in the most massive galaxies. 

• At all redshifts, the covering factor at 2J?vir is < 40% for 
A^Hi > lO^"' cm"^ < 10% forLLSs, and < 1% forDLAs. The Hi 
gas covering factor within 27?vir appears to decrease slowly from 
z ~ 4 to z ~ 1.5, and the covering factor within a fixed aper- 
ture of 10 arcsec (~ 70 — 85 kpc) seems not evolving or slowly 
increasing with redshift. However, the small sample and the non- 
homogeneous sampling of the halo mass function at all the red- 
shifts prevent us from deriving robust scaling relations. In these 
simulations, a certain degree of self-similarity is found in the col- 
umn density probability distribution function for neutral hydrogen 
and the time evolution of the different covering factors appears to 
follow the universal expansion. 

• Cold streams mostly appear as LLSs and contribute signif- 
icantly to the column density distribution function in between 
A^Hi = 10^^ — 10^* cm~'^. Massive disks, as described by these 
high resolution simulations, shape the column density distribution 
function for very high column density DLAs. These models, espe- 
cially if molecules are included, can reproduce the observed knee 
at 10^^'® cm~^ in the f{Nni,X). Conversely, there is a deficiency 
of systems between 10^*^ "' — 10"^" cm^'^, likely due to the missing 
low mass systems in our simulated sample and to finite numerical 
resolution which does not capture a clumpy medium. The simulated 
f{Nni,X) is very shallow below lO^'^ cm~^, indicating that cold 



streams in proximity to massive galaxies are not the only sources 
of opacity for Lyman limit photons in the Universe; an additional 
contribution from the IGM is required. Finally, in agreement with 
observations, the simulated f{NB_i,X) is not evolving in shape 
with redshift, due to a nearly self-similar distribution of neutral hy- 
drogen in halos. 

• Despite the small covering factor, comparing the observed and 
simulated incidence of absorption line systems, it appears that mas- 
sive galaxies and cold streams are responsible for at least 30% of 
the absorbers detected in the foreground of QSOs. In fact, due to a 
large enough nvunber density of galaxies in the vmiverse, the cross 
section required to match the observed incidence is only a factor of 
a few above the mean values in this sample. A modest extrapola- 
tion towards lower halo masses can fully reproduce the observed 
population. To the extent that cold streams are indeed the com- 
mon mode of accreting gas in massive galaxies at high redshift, as 
emerges from the cosmological simulations, these streams have al- 
ready been detected in absorption line surveys, primarily as metal- 
poor LLSs. 

• In our simulations, most of the cross section within 2J?vir is 
occupied by metal poor gas (Z ~ 10~^Zq), that has been partly 
eirriched by previous episodes of star formation. DLAs and SLLSs 
exhibit broad metallicity distributions, with a peak at ^ W~^Zq 
and Z ~ 10'^'^Zq, respectively. Metallicity is a valuable tool 
to disentangle gas that is associated with galaxies or winds (with 
metallicity distribution skewed above Z ~ 10~^^^ZQ) from LLSs 
in the cold streams (with the lowest level of enrichment). 

• Because of the partial covering factor of neutral gas, com- 
posites of Lya absorption lines are not black at the line cen- 
ter, even within 40 kpc from the central galaxy. The typical line 
width outside the central galaxy and its immediate neighborhood is 
~ 200 - 300 krri/s FWHM reflecting the bulk inflow velocity in 
cold streams. Our current simulations reproduce the Lya equiv- 
alent width distribution observed in the circumgalactic medium 
of Lyman break galaxies. Conversely, the predicted optical depth 
of metal lines is suppressed by more than 5 orders of magnitude 
and the predicted equivalent width of SillA1260 is systematically 
lower than the observed values. An increase in the metallicity of the 
streams is not enough to match the observations. Another mech- 
anism is required, such as massive outflows beyond the existing 
winds in the current simulations. 

From the systematic analysis of our simulations and a care- 
ful comparison with observations, we conclude that the theoretical 
predictions of cold gas accretion along filaments of the cosmic web 
are consistent with present day observations. This conclusion is not 
driven by a null claim that the cold streams are undetectable in ab- 
sorption, which would have been rather unsatisfactory. In fact, our 
current simulations, with all the caveats previously discussed, indi- 
cate that cold streams contribute significantly to the incidence and 
kinematics of absorbers observed at redshift z ~ 2 — 3. 

Progress can now be made in two directions. While this work 
offers a first quantitative analysis of detectability of cold streams 
in absorption, simulations representative of a larger spectrum of 
masses will facilitate robust statistical comparisons with observa- 
tions. At the same time, improved modeling of outflows in sim- 
ulations will allow a more detailed exploration of the interaction 
between inflowing and outflowing gas. This is important to bet- 
ter assess the cross section distribution of neutral hydrogen, a key 
quantity for the visibility of cold streams in absorption. 

Observationally, it should be possible to identify candidates 
of cold streams by searching the increasing spectroscopic data sets 
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of quasars and galaxy pairs. Despite the difficulties of measuring 
metallicity in LLSs, one can investigate their kinematics in detail to 
disentangle between metal rich gas that is outflowing at velocities 
> 600 km/s from gas that is inflowing at velocities ~< 200 km/s 
with metallicity < O.OIZq. 

Obser ved metal-poor LLSs are suggestive of the c old mode of 
inflow (e.g. IProchter et allbOlOl: iRibaudo et al]|2011bl) and a con- 
firmation of this phenomenon may arise if a nearby galaxy is ob- 
served to be fed by cold inflows. However, at low redshift the cold 
mode is e xpected to be limited to g alaxies with stellar mass below 
10"Mo jPekel & Bimboiml2006l) , and to accrete from a wide an- 
gle at a lower column density. This is significantly different from 
the cold, narrow, dense streams that are predicted at high redshift 
and only direct evidence for cold inflow at high redshift would be 
a confirmation of these theoretical predictions. In future, provided 
reliable estimates from simulations of observables such as the col- 
umn density distribution, metallicities and kinematics, models of 
cold gas accretion can be also tested statistically against observa- 
tions. 

Following these paths, it may be possible to uncover the 
modes of gas accretion in high redshift galaxies, and their interplay 
with outflows. 
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APPENDIX A: A DETAILED COMPARISON OF THE RT 
CALCULATIONS 

In this appendix, we provide a more detailed comparison of the re- 
sults from the different radiative transfer calculations. At first, we 
discuss the effects that the UVB and local sources have on the col- 
umn density distribution. Then, we quantify the hydrogen density 
threshold above which the gas is self-shielded and we provide a 
numerical approximation useful to correct the CIE approximation 
in the absence of the UVB. Finally, we examine the impact that 
photoionization has on the neutral mass of galaxies and satellites. 



Al Effects on the hydrogen column density distribution 

In Figure lAII we highlight the effect that different ionization 
sources have on the H I column density distribution by showing the 
number of AMR cells in bins of A'hi from an arbitrary projection 
of MW3 at 2 = 2.3. The different RT models have been normal- 
ize with respect to the STAR model (black horizontal solid line). 
The CIE model (blue dash-dotted line) offers an upper limit to the 
neutral fraction at all the column densities. The CIE approximation 
holds only in highly optically thick regions (TVhi > 10^" cm~^), 
while at lower column density overpredicts by a large amount the 
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Figure Al. Compaiison of the neutral hydrogen column density in the 
STAR model (black line), the CIE model (blue dash-dotted line), and the 
UVB model, with and without dust (red dashed line and open squares). The 
STAR model corrected for molecular hydrogen is show with open circles. 
For each RT calculation, we plot the number of projected cells in a given in- 
terval of Afjj I, in compaiison to the STAR model. The CIE approximation 
holds only in highly optically thick regions. At high column densities, neu- 
tral hydrogen is further ionized by local sources and depleted by molecules. 



neutral fraction in comparison to the UVB case (red dashed line). 
Photoionization cannot be neglected, particularly given the fact that 
streams are not entirely self-shielded, as evident from Figure[3] 

Comparing the UVB with the STAR model, it appears that 
radiation from local sources affects primarily the high column den- 
sities (A''hi > 10^^ cm^^) which are completely self-shielded from 
the UVB. Since most of the ionizing photons are absorbed locally, 
the effect of local sources is modest below TVhi ^ 10'^" cm^'^ and 
almost negligible below TVhi ~ 10^* cm~^. Note that quantities 
shown in Figure IaTI are normalized and not proportional to the area 
covered by gas at a given column density. In fact, radiation from 
local sources that leaks form the central disk can affect the H I col- 
umn density distribution function in the interval 10^* — 10^'^ cm~^, 
when weighted by area. Dust produces little or no effect to the UVB 
model (compare the red dashed line with open squares). 

Although not included in the our RT calculation, in Figure IaTI 
we explore how molecules affect the shape of the hydrogen column 
density distribution. When atomic hydrogen is shielded from UV 
radiation, molecules can form, lowering the effective H I column 
density. To investigate the importanc e of this effect, we c o mpute 
the molecular fraction following the iMcKee & Krumholzl j2O10h 
formalism. This analytic model re produces to first order the molec- 
ular gas fracti on in nearbv spirals (|Kruniholz et al .120091) and metal 
poor dwarfs dFumagalli et al.l 201ol) and approximates t he results 
from numerical simulation dKrumholz & GnedinI 1201 ih . In each 
AMR cell, we estimate the molecular gas fraction as a function of 
the total neutral gas column density Egas = 2mpA^H2 + rnpNm 
and metallicity. Since this model is designed to describe the molec- 
ular fraction in individual clouds, we correct for smearing due to 
resolution by assuming Edoud = cEgas with c = 5 for all the 
AMR cells with sizes above 200 pc and c = 1 elsewhere. Although 
not a dominant correction, the inclusion of molecules lowers the H I 
column densities above A^hi ~ 10^^ cm~^ (open circles). For this 
reason, the values of A^h i used in our analysis should be regarded 
as upper limits at the highest density. This column density thresh- 
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Figure A2. Neutral fraction as a function of tlie liydrogen volume density for the CIE model (top left), the UVB model (top right) and the STAR model (bottom 
left) in MW3 at 2 ~ 2.3. The color scale is proportional to the number of AMR cells in each bin, with the most populated region in dark gray. The arithmetic 
mean in bin of log rtj] for the UVB model in individual galaxies is shown in the bottom light panel (gray lines). The central values common to all the galaxies 
is also superimposed (black sohd line). Although with significant scatter, this distribution can be used as a crude approximation to the UVB model. RT with 
local sources is required to capture the intiinsic scatter in the neutral fraction above the ^^H ~ 0.01 — 0.1 cm^'^, where gas becomes shielded from the UVB. 



old is in good agreement with the knee in the observed DLA col- 
umn density distribution JProchaska et al.ll2005l : iNoterdaeme et al.l 
l2009l) that, indeed, is i nterprete d as the transition between atomic 
and molecular gas (e.g. lZwaair& Prochaska.2006) . 

A2 Self shielding and numerical approximations 

Figure lA2l shows the neutral fraction xui as a function of total gas 
volume density nn, for the CIE model (top left panel), the UVB 
model (top right panel), and the STAR model (bottom left panel), 
again for MW3 at z — 2.3. In presence of the UVB, hydrogen is 
highly ionized below nn — 0.01 cm^'^, whereas it is fully neutral 
above nn = 0.1 cm^'^. A similar volume density for self-shielding 
has been consistently found with different numerical te chniques 
(e.g. iRazoumov et al.|[2006l : iFaucher-Giguere et al.ll201ol) . justify- 
ing the use of a fixed threshold above which the UVB can be sup- 
pressed while running the simulations. 

Additional photo-heating, currently neglected in our RT cal- 
culations, will produce a smoother transition between ionized and 
shielded regions due the enhanced collisional ionization in partially 
shielded cells. Also, both the UVB and CIE models do not account 



for the presence of ionizing sources that increase the dispersion of 
xm (bottom left panel in Figure |A2j. A significant fraction of the 
gas above nu — 0.01 — 0.1 cm^'^ remains fully neutral, but local 
fluctuations in the ionization state of up to 5 order of magnitude 
can be found. Indeed, a locus of points with nu > 0.01 cm^'^ and 
xhi < 0.01 is visible in Figure lA2l Since this gas mostly resides in 
the central disk and in satellites where the majority of the star for- 
mation is present, the inclusion of ionizing photons seem to be an 
essential ingredient in models of galaxy formation. Ionizing pho- 
tons from stars are absorbed mostly in the inner part of the galaxy 
by dust and gas and there is only a minor contribution of the ioniza- 
tion below nu = 0.01 cm"'', in qualitative agreement with current 
constraints on the escape fraction of ionizing radiation. However, 
the degree of precision required in comparing models and obser- 
vations requires the inclusion of this RT effects on the neutral gas, 
commonly neglected or treated with some approximations. 

Although below nn = 0.01 cm^'^ the scatter in xhi at a given 
volume density is large, most of the AMR cells accumulate in a well 
defined locus of points. Due to a self-similarity in the probability 
distribution function of the hydrogen column density, there is little 
scatter for different galaxies both with mass and redshift, as seen 
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Figure A3. Comparison of the neutral and total hydrogen mass enclosed within /?vir tor the three RT models. In each panel, main halos are color coded as 
in Figure [T] Satellites are displayed with gray circles. Dotted lines indicate a 1:1 relation, while dashed lines show a robust linear regressions for the main 
halos only. In massive galaxies, only 20-25% of the total hydrogen mass is neutral, with local sources and electron collisions being responsible for most of the 
ionization. The UVB has no effects on the neutral mass due to nearly complete self-shielding at the highest hydrogen volume densities. 



from the bottom right panel of Figure IA2I where we plot the arith- 
metic mean of a^ni for the UVB model (gray hnes for individual 
galaxies). This invariance not only corroborates the assumption of 
a fixed threshold above which the UVB is suppressed, but also al- 
lows one to assume a lower limit on the hydrogen volume density in 
systems that host neutral gas, such as DLAs. Further, we can derive 
a crude approximation to the UVB model. In the CIE model, fully 
neutral or marginally ionized gas is found even below riH = 0.01 
cm~^, with an error up to several order of magnitudes in compari- 
son to the UVB model. Using the fact that the mean xm as a func- 
tion of riH is sensitive to the intensity of the UV background only, 
we derive a relation xm = xm{nYi) that improves the CIE model 
in presence of a UVB. We obtain such relation by combining all the 
arithmetic means computed in individual galaxies with a geometric 
mean (see Figure lA2l and Table |B3t . 

A3 Neutral hydrogen mass 

For each RT model, we measure the total Hi mass enclosed 
in the virial radius of the main halos and satellites (see Table 
IB2I (. A comparison is provided in Figure I A3 1 In the top panel, 
we show the relation between the total and the neutral hydro- 



gen mass in the STAR model. In massive halos above 10 M©, 
less than ~ 1/5 of the hydrogen mass is neutral. A robust lin- 
ear regression in logarithmic space for the main halos only gives 
log A/hi = 0.03 + 0.93 log Mh. Less massive satellites cover a 
wider range in H I mass, from almost neutral to highly ionized. The 
CIE and UVB models predict nearly identical masses above 10^ 
M0 (logMHi.ciE = 0.02 + 1.00 log A/hi, UVB for the main ha- 
los), while in the satellites the CIE model predicts larger neutral 
hydrogen masses. H I masses from the STAR model are systemat- 
ically lower than those in the CIE or the UVB model above 10^ 
M0, up to 50%-60% above 10^° M© (log Afm ,STAR = 0.41 + 
0.93 log A/hi, CIE and log A/hi, star = 0.45 + 0.93 log A/hi,uvb 
for the main halos). Satellites have a larger scatter. 

The explanation for this behavior naturally emerges from the 
previous considerations. In the main halos, most of the mass is in 
high density regions where the gas is self-shielded from the UVB. 
Photoionization from local sources and, to second order, electron 
collisions regulate the amount of mass in neutral hydrogen. For 
satellites, the observed scatter is due to the different SFRs in these 
systems and on the different distances from the central disks. Also, 
at lower masses, systems are more vulnerable to ionization from the 
UVB. This behavior is qualitatively in agreement with the Simula- 
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Figure A4. Bottom panel: column density distribution function for the 
STAR model (black solid line), for a model with increased gas tempera- 
ture (blue dash line) and a model at lower temperature (red dotted line). Top 
panel: the ratio of the distribution for increased (decreased) temperature 
and the STAR model is displayed with a dashed (long dashed) line. Identity 
between the models is marked with a dotted line. Additional heating pro- 
duces higher ionization in the interval log A^ui = 18-3 — 19.8 cm~^, but 
the covering factor of optically thick gas is nearly unperturbed, owing to a 
rising f{Niii,X) towards lower column densities. Lowering the temper- 
ature, little differences (less than a factor of two) are visible since the UV 
radiation is the dominant source of ionization in unshielded regions, while 
shielded regions are already cold and neutral. 

tions bv lGnediiJlioiol) who finds that the Lyman- Werner radiation 
from the UVB is neghgible compared to the local radiation field 
and that the ionizing radiation from the UVB plays a significant 
role only up to the edge of the Hi disks. This result has implica- 
tions for understanding the processes that shape the massive end of 
the H I mass function and stresses further the necessity to account 
for l ocal photoionization as an important feedback mechanism (e.g. 
ICantalupa.20ia) . 



A4 Uncertainties and future improvements 

Although our model significantly improves upon simpler RT cal- 
culations, few approximations contribute the final error budget. For 
example, these simulations do not include the presence of an active 
galactic nucleus and additional ionization from this harder spec- 
trum is neglected in this work. Similarly, we include only pho- 
toionization at 912A due to stars and we do not considering harder 
SEDs that could enhance the ionization via heating. Similarly, He 
recombination emission is neglected. Finally, we ignore the pos- 
sibility that these galaxies could l ie in the proximity of other ion- 
izing sources, such as QSOs (e.g. ICantalupo et al.ll2005h . We note 
the these effects would reduce further the neutral gas fraction and 
enhance the effects of ionizing radiation on the neutral volume and 
column density. 

In the previous sections, we have clearly shown that the in- 
clusion of both the UVB and local sources shapes the distribution 
of gas volume and column density at the low and high densities. 
However, the exact amplitude of these effects is proportional to the 
intensity of the input radiation field which is subject to uncertain- 
ties, as discussed in Section \3A\ Heating from the UVB and lo- 
cal sources is already included in these simulations, but during the 



post-pr ocessing we do not update th e gas temperature (see a discus- 
sion in lCantalupo & Porcianill2oTll) . And a sensible determination 
of the temperature has o bvious im plications for collision ionization 
and emission properties dFaucher- Giguere et al. 2010). 

We try to quantify the impact that an increase in tempera- 
ture has on our results by performing a simple test with MW3 
at 2: ~ 2.3, using the STAR model. After having increased by 
50% the temperature in photo-ionized cells with xm < 0.8, we 
recompute the neutral fraction for CIE (a;Hi,cic). Then, we de- 
rive the projected H l distribution function using an updated neu- 
tral fraction a^ni.now = min(^a;Hi,oida;Hi,cio; a;Hi,oid; a;Hi,cic), de- 
fined in this way to favor high ionization. Results are summa- 
rized in Figure IA4I In the bottom panel, we show the column 
density distribution function for the original STAR model (black 
solid line) and for the model at higher temperature (blue dash 
line). In the top panel we display the ratio of the two distribu- 
tions. Additional heating produces higher ionization in the inter- 
val logiVni = 18.3 - 19.8 cm'^, yielding a flatter /(7Vhi,X). 
The cross section of SLLSs is reduced by a factor of ~ 2, in- 
creasing even further the discrepancy with observations. The DLA 
cross section is decreased by no more than 20% and the effect on 
the cumulative covering factor of LLSs is small, owing to a ris- 
ing f(Nni,X) towards lower column densities. This simple test 
suggests that a moderate increase is temperature does not have a 
major impact on the result presented in this paper. Conversely, if 
cooling is more effective than what assumed in these simulations, 
hydrogen is less ionized by electron collisions. In Figure lA4l we 
shown the hydrogen distribution (red dotted line) for a model in 
which the temperature is reduced to 50% of its value in all the cells 
with TiH > 0.05 cm~^, i.e. partially or totally shielded. Gas be- 
low this density is highly photoionized and the ionization state is 
insensitive to the gas temperature. In this case, the ionization frac- 
tion is assumed to be a^ni.now = a;Hi,oid for nu < 0.05 cm~^ 
and a;Hi,ncw = a^ni.cic for nn ^ 0.05 cm~^. Since gas above 
riH ~ 0.05 cm~^ is already at low temperature, little difference 
(less than a factor of two) is visible in the final H I distribution. 
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Table Bl. Summary of the baryonic and dark matter properties of the 7 galaxies used in this study. 



Redshift 


MWl 


MW2 


MW3 


MW4 MW5 


MW8 


MW9 


MWl 


MW2 


MW3 


MW4 MW5 


MW8 


MW9 








Virial Mass (lO" Mq) 










Virial Radius (kpc) 






4.00 


0.11 


0.91 


0.29 


1.26 6.09 


0.62 


0.19 


15.00 


30.00 


20.50 


33.50 56.50 


26.25 


18.25 


3.55 


0.26 


1.51 


0.37 


1.73 9.58 


0.73 


0.33 


21.75 


39.00 


24.50 


41.00 72.25 


30.50 


23.50 


3.17 


0.53 


3.14 


1.06 


2.73 11.70 


1.13 


0.48 


29.75 


57.00 


38.00 


51.75 84.25 


38.50 


29.00 


2.85 


1.13 


4.19 


1.62 


4.23 


1.52 


1.22 


41.75 


64.50 


47.25 


64.50 


46.00 


42.75 


2.57 


1.68 


5.04 


2.66 


5.81 


1.80 


1.31 


51.50 


74.25 


60.50 


77.25 


52.25 


47.00 


2.33 


3.34 


5.92 


3.55 


9.21 


2.31 


1.29 


68.75 


83.25 


70.50 


96.50 


60.75 


50.00 


2.13 


6.49 


7.64 


5.73 


12.01 


2.54 


1.32 


91.25 


97.25 


88.00 


1 12.25 


66.75 


53.50 


1.94 


8.83 


9.12 


8.06 


15.17 


2.79 


1.75 


106.75 


108.75 


104.00 


129.00 


73.00 


62.50 


1.78 


10.20 




9.27 




2.85 


2.87 


118.50 


_ 


114.75 




77.75 


77.50 


1.63 


11.47 




11.81 




2.85 


3.00 


129.50 


_ 


131.25 




81.75 


82.75 


1.50 


13.03 




12.59 




2.86 


3.39 


141.75 


_ 


140.25 




85.50 


90.50 


1.38 


13.69 




13.52 




5.10 


4.31 


150.50 


_ 


150.00 




108.00 


102.25 






Star Formation Rate (Mq yr ^) 










Stellar Mass (10^° Mq) 






4.00 


0.94 


37.79 


3.74 


21.19 49.49 


4.23 


2.02 


0.03 


1.62 


0.20 


1.22 6.25 


0.63 


0.13 


3.55 


4.61 


58.25 


2.70 


28.66 140.61 


5.70 


2.88 


0.17 


3.12 


0.26 


1.69 10.66 


0.80 


0.28 


3.17 


4.42 


148.23 


14.44 


15.77 76.14 


4.48 


3.12 


0.31 


7.52 


0.86 


2.79 13.32 


1.31 


0.41 


2.85 


14.76 


157.48 


5.54 


38.11 


40.65 


6.98 


0.87 


11.16 


1.14 


3.91 


1.70 


1.27 


2.57 


9.65 


104.36 


25.95 


48.22 


4.73 


8.28 


1.14 


14.24 


1.76 


4.94 


1.85 


1.48 


2.33 


35.71 


99.77 


30.49 


44.21 


7.16 


7.14 


2.25 


17.51 


2.63 


8.69 


2.21 


1.61 


2.13 


52.47 


143.23 


46.77 


73.75 


5.03 


2.45 


5.54 


22.26 


4.34 


12.07 


2.42 


1.70 


1.94 


62.03 


187.86 


71.46 


77.22 


13.63 


2.38 


8.19 


27.48 


7.14 


15.38 


2.96 


1.94 


1.78 


108.52 




84.00 




6.26 


5.44 


10.57 


_ 


8.90 




3.04 


3.19 


1.63 


20.71 




50.19 




3.42 


17.68 


11.43 


_ 


12.10 




3.09 


3.46 


1.50 


25.65 




39.10 




3.61 


4.25 


12.80 


_ 


13.15 




3.20 


3.98 


1.38 


23.19 




45.64 




5.31 


6.77 


14.19 


_ 


14.38 




6.42 


5.49 








Dark Matter Mass (lO" Mq) 










Hydrog 


en Mass (10^" Mq) 






4.00 


0.10 


0.66 


0.25 


1.07 5.19 


0.53 


0.16 


0.080 


0.879 


0.132 


0.663 2.761 


0.291 


0.140 


3.55 


0.23 


1.06 


0.33 


1.49 8.13 


0.61 


0.29 


0.167 


1.431 


0.174 


0.791 3.921 


0.361 


0.182 


3.17 


0.46 


2.04 


0.93 


2.32 9.86 


0.95 


0.41 


0.320 


3.517 


0.480 


1.265 5.041 


0.566 


0.232 


2.85 


0.98 


2.66 


1.42 


3.62 


1.29 


1.05 


0.619 


4.138 


0.906 


2.156 


0.635 


0.494 


2.57 


1.46 


3.24 


2.30 


5.01 


1.55 


1.11 


0.996 


3.803 


1.823 


3.057 


0.614 


0.447 


2.33 


2.88 


3.72 


3.07 


7.84 


2.01 


1.09 


2.390 


4.502 


2.249 


5.065 


0.862 


0.419 


2.13 


5.54 


4.65 


4.95 


10.23 


2.20 


1.11 


3.990 


7.622 


3.453 


5.778 


1.058 


0.437 


1.94 


7.54 


5.43 


6.92 


12.88 


2.39 


1.48 


4.648 


9.399 


4.181 


7.590 


1.109 


0.767 


1.78 


8.70 




7.98 




2.43 


2.42 


4.449 




4.009 




1.136 


1.312 


1.63 


9.85 




10.11 




2.42 


2.53 


4.791 




4.926 




1.162 


1.255 


1.50 


11.17 




10.77 




2.42 


2.86 


5.763 




5.097 




1.177 


1.382 


1.38 


11.68 




11.54 




4.27 


3.59 


5.920 




5.413 




1.895 


1.701 



Table B2. Summai y of the neutral gas properties of the 7 galaxies used in this study. 



Redshift 


MWl 


MW2 


MW3 


MW4 


MW5 


MW8 


MW9 


MWl 


MW2 


MW3 


MW4 


MW5 


MW8 


MW9 








HI Mass (IQl'' M© 


)-UVB 










HI Mass (10"' Mq 


) - STAR 






4.00 


0.036 


0.565 


0.043 


0.326 


1.124 


0.158 


0.068 


0.026 


0.322 


0.026 


0.215 


0.637 


0.158 


0.059 


3.55 


0.072 


0.926 


0.054 


0.271 


1.423 


0.186 


0.079 


0.041 


0.464 


0.023 


0.124 


0.891 


0.145 


0.060 


3.17 


0.103 


2.210 


0.110 


0.389 


1.147 


0.302 


0.081 


0.069 


0.940 


0.049 


0.270 


0.589 


0.252 


0.057 


2.85 


0.190 


2.631 


0.295 


0.676 




0.218 


0.187 


0.108 


1.281 


0.217 


0.410 




0.091 


0.187 


2.57 


0.285 


2.085 


0.751 


1.071 




0.194 


0.120 


0.176 


1.169 


0.434 


0.651 




0.153 


0.078 


2.33 


0.805 


2.230 


0.869 


1.963 




0.277 


0.074 


0.490 


1.081 


0.522 


1.494 




0.216 


0.044 


2.12 


1.454 


3.950 


1.195 


1.548 




0.379 


0.115 


0.901 


1.983 


0.763 


0.977 




0.379 


0.087 


1.94 


1.537 


4.902 


1.193 


1.902 




0.372 


0.230 


0.850 


2.427 


0.498 


1.252 




0.248 


0.230 


1.78 


0.794 




0.922 






0.476 


0.488 


0.318 




0.323 






0.379 


0.415 


1.63 


0.745 




0.885 






0.544 


0.452 


0.441 




0.369 






0.478 


0.324 


1.50 


1.207 




1.048 






0.579 


0.385 


0.810 




0.591 






0.487 


0.382 


1.38 


1.202 




0.818 






0.922 


0.560 


0.839 




0.407 






0.798 


0.559 
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Table B3. Geometric mean of the Ihi vs njj relation from the UVB model. 



lognH(cm ^) 


log XHI 


lognH(cm ^) 


log XHI 


lognH(cm ^) 


log XHI 


lognnCcm ^) 


log XHI 


-4.750 


-5 51 2F4-00 


-2.950 




-1.150 


-2 522F-ni 


0.650 


-5 050F-03 


-4.650 




-2.850 




-1 050 




750 




-4.550 


-5 496F+0n 


-2.750 


-3 506F+0n 


-0.950 


-1.192E-01 


0.850 


-3 391F-03 


-4.450 


-S 475F+0fl 


-2.650 


-3 279F+0n 


-0.850 




0.950 


-2 770F-03 


-4.350 


-5 445F+nn 


-2.550 


-3 n54F+on 


-0.750 


-8 fi07F-02 


1.050 


-7 3n7F-03 


-4.250 


-5.448E-I-00 


-2.450 


-2.820E+00 


-0.650 


-7.321E-02 


1.150 


-1.850E-03 


-4.150 


-5.401E-I-00 


-2.350 


-2.579E+00 


-0.550 


-6.224E-02 


1.250 


-1.530E-03 


-4.050 


-5.360E-I-00 


-2.250 


-2.323E+00 


-0.450 


-5.388E-02 


1.350 


-1.256E-03 


-3.950 


-5.293E-I-00 


-2.150 


-2.050E+00 


-0.350 


-4.524E-02 


1.450 


-9.829E-04 


-3.850 


-5.212E-I-00 


-2.050 


-1.764E+00 


-0.250 


-3.801E-02 


1.550 


-6.561E-04 


-3.750 


-5.128E-I-00 


-1.950 


-1.432E+00 


-0.150 


-3.123E-02 


1.650 


-4.945E-04 


-3.650 


-5.020E-I-00 


-1.850 


-1.098E+00 


-0.050 


-2.520E-02 


1.750 


-3.221E-04 


-3.550 


-4.906E-I-00 


-1.750 


-8.472E-01 


0.050 


-2.037E-02 


1.850 


-2.372E-04 


-3.450 


-4.775E+00 


-1.650 


-6.688E-01 


0.150 


-1.579E-02 


1.950 


-1.611E-04 


-3.350 


-4.631E+00 


-1.550 


-5.363E-01 


0.250 


-1.247E-02 


2.050 


-1.826E-04 


-3.250 


-4.479E+00 


-1.450 


-4.328E-01 


0.350 


-1.015E-02 


2.150 


-8.635E-05 


-3.150 


-4.321E+00 


-1.350 


-3.546E-01 


0.450 


-7.837E-03 






-3.050 


-4.155E-I-00 


-1.250 


-2.978E-01 


0.550 


-6.192E-03 







Table B4. DLAs and LLSs covering factor from the STAR model. The top table lists values for the streams and galaxies, while the bottom one lists values for 
the streams alone. 



Redshift 


MWl 


MW2 


MW3 


MW4 MW5 


MW8 


MW9 


MWl 


MW2 


MW3 


MW4 MW5 MW8 


MW9 




DLA covering 


factor in 


% at 2i?vir (streams and galaxies) 


LLS covering factor in 


% at 2i?vir (streams and galaxies) 


4.00 


1.78 


3.47 


1.09 


2.68 1.42 


1.57 


1.55 


10.80 


12.54 


6.70 


12.19 12.43 11.20 


11.72 


3.55 


1.03 


3.65 


0.79 


1.35 1.22 


0.85 


1.02 


7.92 


9.02 


6.24 


8.56 11.11 9.35 


8.27 


3.17 


0.94 


1.95 


0.60 


1.76 1.08 


0.88 


0.76 


8.03 


5.08 


2.35 


13.12 8.68 4.71 


6.41 


2.85 


0.79 


2.65 


1.42 


1.86 


0.28 


0.62 


4.88 


5.81 


10.41 


11.03 - 1.69 


5.36 


2.57 


1.12 


2.84 


0.88 


1.60 


0.31 


0.29 


10.25 


7.13 


6.52 


10.30 - 4.08 


2.98 


2.33 


1.37 


2.22 


1.19 


1.54 


0.31 


0.49 


9.87 


6.89 


8.92 


9.11 - 3.80 


4.72 


2.12 


1.07 


2.23 


0.75 


1.00 


0.59 


0.82 


6.64 


6.79 


5.74 


7.03 - 4.03 


5.06 


1.94 


0.62 


2.07 


0.38 


0.72 


0.28 


0.73 


4.91 


6.85 


3.19 


4.85 - 2.03 


5.09 


1.78 


0.23 




0.17 




0.25 


0.66 


1.98 




1.83 


1.63 


4.20 


1.63 


0.36 




0.21 




0.51 


0.30 


4.25 




2.84 


2.46 


2.59 


1.50 


0.35 




0.27 




0.51 


0.35 


2.97 




2.82 


2.87 


5.04 


1.38 


0.26 




0.32 




0.48 


0.44 


1.84 




2.40 


1.66 


4.21 






DLA covering factor in % at 2iJvir (streams only) 






LLS coverin 


g factor in % at 2/?vir (streams only) 




4.00 


1.47 


2.94 


0.82 


2.41 0.77 


0.95 


1.29 


10.55 


12.05 


6.30 


11.63 11.82 10.60 


11.42 


3.55 


0.70 


3.02 


0.52 


0.98 0.68 


0.34 


0.74 


7.46 


8.28 


5.78 


7.95 10.41 8.75 


7.83 


3.17 


0.62 


1.48 


0.36 


1.28 0.39 


0.33 


0.57 


7.51 


4.41 


1.87 


12.45 7.85 3.96 


6.09 


2.85 


0.41 


2.24 


1.00 


1.27 


0.05 


0.37 


4.17 


5.23 


9.67 


10.24 - 1.13 


4.85 


2.57 


0.69 


2.40 


0.32 


0.91 


0.03 


0.15 


9.73 


6.58 


5.51 


9.34 - 3.40 


2.49 


2.33 


0.70 


1.89 


0.54 


0.95 


0.08 


0.15 


8.99 


6.42 


7.88 


8.29 - 3.13 


4.04 


2.12 


0.50 


1.83 


0.32 


0.50 


0.15 


0.26 


5.66 


6.31 


5.10 


6.26 - 3.31 


4.21 


1.94 


0.22 


1.75 


0.16 


0.44 


0.03 


0.37 


4.23 


6.41 


2.60 


4.38 - 1.56 


4.50 


1.78 


0.07 




0.04 




0.03 


0.34 


1.62 




1.44 


1.02 


3.42 


1.63 


0.08 




0.07 




0.24 


0.05 


3.64 




2.47 


1.98 


2.04 


1.50 


0.07 




0.07 




0.15 


0.12 


2.30 




2.31 


2.22 


4.30 


1.38 


0.02 




0.14 




0.28 


0.24 


1.34 




1.97 


1.19 


3.64 
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Table B5. /{Nm, X) from the STAR model witin -Rvir, 2_Rvir and in the streams alone. 



logATni (cm-2) 


Rvii 




2i?vir 




Streams 




13.5 


2.22 X 10" 


-15 


8.18 X IQ- 


-15 


8.30 X 10- 


-15 


14.5 


2.61 X 10" 


-16 


1.73 X 10" 


-15 


1.74 X 10" 


-15 


15.5 


2.16 X 10" 


-17 


1.35 X 10- 


-16 


1.36 X 10- 


■16 


16.5 


1.94 X 10" 


-18 


5.98 X 10" 


-18 


6.06 X 10- 


-18 


17.5 


9.10 X 10- 


-20 


3.03 X 10- 


-19 


3.00 X 10- 


-19 


18.5 


2.91 X 10" 


-21 


9.72 X 10" 


-21 


9.52 X 10- 


■21 


19.5 


2.64 X 10- 


-22 


7.70 X 10- 


-22 


7.58 X 10- 


■22 


20.5 


2.65 X 10" 


-23 


6.43 X 10- 


-23 


5.46 X 10- 


-23 


21.5 


1.38 X 10" 


-24 


2.55 X 10- 


-24 


1.17 X 10- 


-24 


22.5 


7.80 X 10" 


-27 


1.26 X 10- 


-26 


1.78 X 10- 


-28 


23.5 


9.92 X 10" 


-30 


1.24 X 10- 


-29 






24.5 


4.08 X 10" 


32 


4.08 X 10" 


-32 







Table B6. Average DLA cross sections within i?,vir for aU the halos at all redshifts. 



Halo Mass 


Total cross section 


Central galaxy only 


(M©) 




(kpc2) 




(kpc2) 




Mean 


Standard deviation 


Mean 


Standard deviation 


3.6 X IQio 


46 


24 


30 


21 


8.6 X 10^° 


138 


93 


82 


58 


1.6 X 10" 


209 


170 


104 


86 


2.6 X 10" 


312 


225 


182 


145 


3.9 X 10" 


723 


396 


251 


184 


6.8 X 10" 


623 


498 


215 


252 


1.1 X 10^2 


522 


252 


216 


163 
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Table B7. Equivalent widths for stacked Lya and Si IIA1260 absoiption profiles at different impact parameters. 
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